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INTRODUCTION 

The  advent  of  the  space  age,  and  the  proposals  for  nuclear  rockets  which 
have  naturally  followed,  has  presented  perplexing  physical  situations  which 
have  not  been  previously  considered.   One  of  the  more  prominent  of  these  is 
that  of  the  heating  of  materials  by  the  attenuation  of  a  beam  of  gamma  photons 
of  high  flux  intensity. 

Such  a  physical  situation  would  be  found  in  the  structural  members  of  a 
nuclear  rocket  which  are  in  the  near  proximity  of  the  reactor.   Although  the 
problem  of  non  uniform  heat  generation  is  somewhat  unique,  it  is  really  the 
unusual  boundary  conditions  which  give  rise  to  this  new  series  of  problems  which 
are  dissimilar  from  those  which  have  been  conventionally  considered.   Conse- 
quently, the  mathematical  models  needed  to  analyze  such  systems  have  also  need- 
ed alteration. 

Two  properties  of  special  interest  in  such  systems  are  the  temperature 
and  the  thermal  stresses  in  the  materials  subjected  to  this  high  intensity 
gamma  photon  flux.   Structural  members  obviously  will  have  their  utility  seri- 
ously impaired  if  the  temperature  of  the  material  is  allowed  to  approach  the 
melting  point  or  if  the  stresses  arising  in  the  material  rise  to  such  propor- 
tions .that  the  member  tends  to  buckle. 

In  keeping  with  the  physical  conditions  under  consideration,  the  avail- 
ability of  certain  modes  of  heat  transfer  can  be  considered  only  in  specific 
instances.  The  most  prominent  evidence  of  this  is  the  absence  of  convection 
in  regions  of  extremely  low  atmospheric  density  such  as  would  occur  in  outer 
space.  Consequently,  the  primary  manner  in  which  heat  may  be  carried  from  a 
body  is  first  through  conduction  to  the  surface  and  then  by  thermal  radiation 
from  the  surface  to  the  surroundings. 


Since  a  limited  amount  of  information  is  available  in  the  literature 
with  regard  to  such  systems,  it  is  necessary  to  construct  mathematical  models 
which  closely  approximate  these  conditions  and  then  attempt  to  find  solutions 
to  the  resulting  equations  and  boundary  conditions  which  are  proposed. * 

A  first  approximation  to  a  feasible  physical  system  is  that  of  an  in- 
finite plate  with  insulation  at  one  face  and  radiative  heat  loss  at  the  other 
face.  Such  a  system  easily  yields  a  solution,  •  especially  in  the  steady  state 
case,  provided  the  flux  is  of  a  constant  intensity.  However,  since  radiative 
heat  loss  presents  a  non  linear  system  of  equations,  it  is  necessary  to  -ap- 
proximate the  true  solution  by  forming  an  equivalent  gray  body  convective 
radiation  coefficient  so  that  a  linear  system  is  created.  Equating  this  to 
the  true  thermal  radiation  heat  loss  determines  the  value  of  the  coefficient 
and  subsequently  the  exact  solution  to  the  resulting  system  of  equations. 

During  the  operation  of  the  reactor,  the  gamma  flux  can  be  considered 
to  be  of  constant  intensity.  However,  after  shutdown  the  flux  becomes  a 
function  of  time  and  the  problem  of  expressing  the  temperature  of  the  plate 
with  regard  to  its  time  dependency  becomes  considerably  more  complex.  This 
complexity  is  significantly  reduced  by  resorting  to  numerical  techniques 
coupled  with  the  use  of  a  high  speed  digital  computer.  The  primary  signifi- 
cance of  the  considerations  involved  in  such  a  situation  is  the  knowledge 
of  the  time  evolved  in  attaining  the  maximum  temperature  of  the  plate  after 
the  shutdown  of  the  reactor  and  also  the  magnitude  of  this  temperature.  Of 
course  the  period  of  operation  of  the  reactor  and  the  type  of  material  being 
considered  are  determining- factors  in  such  considerations. 

In  an  attempt  to  keep  the  temperature  and  the  maximum  thermal  stress 
well  within  the  probable  destructive  limitations,  it  is  possible  to  construct 


a  system  whereby  additional  heat  transfer  can  be  made  available.  The  con- 
struction considered  here  is  that  of  two  parallel  plates  which  would  have  the 
same  total  thickness  as  a  similar  plate  which  has  been  previously  considered . 
The  inside  face  of  the  first  plate  is  insulated  and  the  outside  face  of  the 
second  plate  loses  heat  to  its  surroundings  only  by  thermal  radiation.  Be- 
tween the  two  plates  heat  is  lost  by  convection  such  as  would  occur  if  a  gas 
were  kept  flowing  through  this  channel. 

An  accurate  evaluation  of  this  system  would  involve  an  extremely  detail- 
ed and  complex  mathematical  model.  However,  a  good  indication  of  the  behavior 
of  the  system  can  be  obtained  by  assuming  a  single  convective  coefficient  to 
be  applicable  over  the  entire  face  of  the  plate.  Then  by  varying  this  con- 
vective coefficient  over  a  fairly  wide  range,  the  temperature  response  of  the 
system  can  be  investigated.  The  minimization  of  the  maximum  thermal  stress 
for  each  of  the  two  plates  is  obtained  by  altering  the  size  of  both  plates 
while  keeping  the  total  thickness  constant  until  the  average  temperature  in 
eaoh  of  the  plates  is  equivalent.  This  choice  is  determined  by  the  faot  that 
the  maximum  thermal  stress  is  directly  proportional  to  the  temperature. 


NOMENCLATURE 

A  Area  (ft2) 

c  heat,  capacity  (Btu/lb,°F) 

E  modulus  of  elasticity  (psi) 

h-j ,hg  convection  coefficient  at  the  surface  of  plate  one  and  two 
respectively  in  the  double  plate  iteration  (Btu/hr,ft2,°F) 

hr  equivalent  gray  body  convective  radiation  coefficient  (Btu/hr, f t2,°F) 

H  volumetric  heat  generation  rate  (Btu/hr,  ft-*) 

I  gamma  photon  flux  intensity  (Btu/hr, ft  ) 

k  thermal  conductivity  (Btu/hr, ft, °F) 

L  thickness  of  the  plate  (ft) 

«(.(  )     Laplace  transform  of  the  variable  in  parenthesis 

q/A  heat  flux  (Btu/hr, ft2) 

8  Laplacian  variable 

t  time  (hr;  seconds  in  time  dependent  flux  equation) 

t_  time  of  reactor  operation  (sec) 

s 

At  time  increment  between  two  successive  calculations  (hr) 

T  temperature  (°R) 

V  volume  (ft^) 

x  distance  through  the  plate  (ft) 

Ax  spatial  increment  in  finite  difference  calculations  (ft) 
(=L/nuraber  of  increments) 

A*  coefficient  of  thermal  expansion 

€.  emissivity 

U  temperature  difference  between  the  body  and  the  heat  sink  (°F) 

O  Laplace  transform  of  8 


vj  temperature  difference  between  the  body  and  the  heat  sink  at  the 

nth  time  increment  and  the  Jth  mesh  point  across  the  plate  in  the 
finite  difference  calculations  (°F) 

0{  thermal  diffusivity  (ft2/hr) 

A»       eigenvalue  which  satisfies  a  transcendental  equation  for  a 
particular  set  of  boundary  conditions  (ft"') 

U  linear  energy  absorption  coefficient  (ft"' ) 

y  Poisson's  ratio  (ft /ft) 

O  density  (lb/ft3) 

C  Steffan  -  Boltzmann  constant  (Btu/hr,ft2,°R ) 

(T  (subscripted)  -   thermal   stress   (psi) 

V  real  constant  used  in  the  finite  difference  calculations 

V  Laplacian  operator 


THEORY 

In  order  to  find  the  thermal  stresses  in  any  body,  it  is  necessary  fo 
investigate  the  resultant  temperature  distribution,  preferentially  determined 
as  the  temperature  increase  above  some  base  temperature.   These  results  ar-; 
directly  related  to  the  type  of  heating  to  which  the  body  is  exposed  and  to 
the  boundary  conditions  to  which  it  is  subjected.   The  following  discussion 
considers  a  gamma  photon  heat  source  under  two  basic  conditions: 

1)  A  non  time  dependent  gamma  source,  such  as  would  occur  during  the 
steady  state  operation  of  a  reactor. 

2)  A  decaying  gamma  source,  such  as  would  occur  after  the  shutdown  of 
a  reactor. 

The  partial  differen+ial  equation  describing  the  temperature  distribution 
in  a  body  with  an  internal  heat  source  is* 

kv*T    +   H    =/c'ff  (1) 

where  the  thermal  conductivity,  k,  and  the  density,^  ,  have  been  assumed  as 
being  constant  at  the  average  temperature  of  tho  body.  A  more  convenient  way 


&   •+«■     -4&  (t) 


of  writing  this  for  an  infinite  plate  is 

The  first  case  considered  will  be  that  of  a  steady  state  gamma  source, 
thus  implying  a  heat  source  which  has  only  spatial  dependence.      The  assumption 
is  made  that   the  flat  plate,  which  is  under  consideration,    is  exposed  to  a 
plane  colliraated  beam  of  gamma  photons.      This  is  a  reasonable  assumption  if 


*see  Nomenclature,    p.   4. 


the  contribution  from  capture  gammas  can  be  ignored.   A  few  elementary  cal- 
culations, using  Grotenhuis  (9),  Table  5-A.  show  that  for  a  plate  of  one  inch 
thickness  this  assumption  is  justified,  since  the  capture  gamma  photon  con- 
tribution is  probably  less  than  five  percent  of  +he  total  flux.   The  heat 
generation  term  is  then  written  as 

H  =/<Ie^  (3) 

where  the  absorption  coefficient  has  been  adjusted  to  consider  the  minimal 
effect  of  buildup  (3)1  O^).   The  equation  t0  be  solved  thus  becomes 

with  the  exact  solution  depending  on  1he  boundary  conditions. 

The  first  investigation  is  the  idealized  situation  in  which  the  side 
of  the  plate  exposed  to  the  gamria  photon  flux  is  perfectly  insulated  and  the 
other  side  of  the  plate  loses  heaf  only  by  thermal  radiation.   It.  is  also 
assumed  that  there  is  an  initially  flat  ' emporature  distribution  at  +  he  same 
temperature  as  the  heat  sink  to  which  the  plate  is  radiating.  Applying  these 
conditions  appropriately,  it  is  found  that  the  temperature  distribution  is 
given  by  (see  Appendix  A) 


(5) 


e-*tt 


where  A„  are  the  n  positive  roots  of  the  transcendental  equation 


and  h  is  an  equivalent  gray  body  convec+ive  radiation  coefficient,  such  that 

Timoshenko  and  Goodier  (21 )  show  that  the  thermal  stresses  in  lhe  jr  and 
z   directions  of  an  unrestrained  flat  plate  with  a  temperature  distribution 
along  the  x  axis  are  given  by 


oj  =  <rB  =  -*^>  +  ^p^^UETW^ 


M6 

*0-/z)4(l-7>) 


Ml  (8) 


-A&; — Z5r\  -Xrt-ETC^)^ 


'-Vt 

where  T  represents  the  temporal ure  rise  above  some  base  temperature  at  which 
the  plate  is  in  an  unstressed  condition.   If  this  base  temperature  is  taken 
to  be  the  same  as  the  initial  temperature  of  the  plate,  Eq.(8)  may  be  written 
in  terms  of  8f  derived  above,  Drovided  the  axis  is  shifted  so  as  to  permit 
integration  from  0  to  L.   This  is  easily  accomplished  by  replacing  x  in 
Bq.(8),  by  the  quantity  (x  _  L/2).   The  equation  then  becomes 


<5  = 


-  or  =  ^M-  [-8  +  f  J'  B'4*  +  vz(^'V^-^)flaxj(9) 


Performing  the  indicated  operations  for  the  case  of  the  temperature  dis- 
tribution shown  in  Bq.(5)  yields  (see  Appendix  A) 

m  r* 


Xh*.  »A'Jf  Kx'.Ww  -jfrWWx.L 


00 


where  G^  is  defined  as 

^""  AX-U^/H  KL+Ji'X^L+Mt 

After  shutdown,  the  gamma  flux  becomes  time  dependent.   Lottes  (12) 
gives  this  dependency  as 

it  /    ***      t      tr*       -(^      .fefeso) 

I(t)xI(o)jp.l(^IO]  -0.067(t+*Mo)-o.OOfc5e   -o.oo\^€ 

(12) 

-fo.ifr*t.*iof-o.oa7(t*^ 

where  t  is  the  time  of  reactor  operation  and  t_  is  the  time  after  shutdown, 
all  time  units  being  in  seconds. 

Because  of  this  time  dependency,  the  most  convenient  analytical  solution 
to  this  problem  is  through  the  use  of  the  Laplace  transformation.  However, 
this  leads  to  an  equation  in  the  Laplace  domain  which  is  quite  unwieldy. 
Consequently,  it  was  deemed  advisable  to  resort  to  the  use  of  numerical  tech- 
niques.  Though  such  a  system  is  not  precisely  accurate,  it  can  be  made  to 
give  answers  to  accuracies  as  close  as  desired  by  choosing  variables  of  prop- 
er magnitude.   Expressing  Eq.(4)  in  finite  difference  form,  the  temperature 
at  the  (n  +  l)th  time  increment  at  distance  x  =  jax  from  the  side  of  the 
plate  is  represented  by 


X\     At      / 


(13) 


where  <+»  is  a  real  constant  0  ^  f  <  1 .     If  £  *<+>i  1,    the  system  is  stable 
and  can  be  solved  for  any  choice  of  Ax  and   ^t.     If  0  -  ^f*  -  i»   the  system 


10 


is  stable  only  when 

Z^^t   -£  _JL_  (1*0 

Richtmyer  (18)  shows  that  suoh  a  system  converges  to  the  analytical  solution 
most  rapidly  and  is  always  stable  when 

This  system  yields  a  set  of  equations  which  may  be  put  in  matrix  form  and 
easily  solved  for  the  temperature  distribution  with  the  aid  of  a  high  speed 
digital  computer.  The  system  of  equations  is  described  in  detail  in  Appendix 
B. 

The  final  area  to  be  considered  is  that  of  two  parallel  plates  with  a 
space  between  them.  This  space  is  assumed  to  act  as  a  channel  through  which 
a  flow  of  gas  is  kept  moving.  As  a  first  approximation,  an  element  of  the 
two  plates  is  considered,  between  which  there  is  a  constant  oonveotive  heat 
transfer  coefficient.  The  object  of  such  an  arrangement  would  be  to  keep 
the  total  thiokness  of  the  two  plates  constant  and  then  vary  the  individual 
thickness  in  such  a  way  as  to  equalize  the  average  steady  state  temperatures. 
The  inner  plate  would  be  subjected  to  the  same  boundary  conditions  as  existed 
in  the  oase  of  the  single  plate,  exoept  that  the  radiation  ooeffioient  would 
be  replaced  by  a  true  oonveotion  coefficient.  Examination  of  Bq.(5)  shows 
that  the  resulting  equation  is 

'  (16) 


11 


where  the  subscript  1  refers  to  the  first  plate.  The  second  plate  loses  ueat 
by  convection  from  its  inner  face  and  by  radiation  from  its  outer  face.  The 
resulting  temperature  distribution  in  this  plate  is  found  to  be 

**"-  ft*  vt  -j^  ^Mfcl  +^fc-  %f 


where 


"\ 


V -«.  ft  ♦  u)  *  i#"tfcHt) 

*nd    AsJLs  given  by 

+      \    \     -*.  j^iAa-khz ""  lis) 

and  l£»  the  flux  incident  on  the  second  plate  is  given  by 


> 


(18) 


(19) 


(20) 


12 


Each  of  the  constants  in  Eq.(l6)   and   (17)  have  been  subscripted   to  denote 
which  plate  is  being  analyzed.     Since  each  plate  is  at  a  different  average 
temperature,    the  physical  constants  are  different  for  each  plate. 
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CALCULATIONS 


The  Single  Plate  with  Constant  Gamma  Photon  Flux  and 
Radiative  Heat  Loss 


The  temperature  rise  within  a  flat  plate  is  a  function  of  the  energy- 
absorbed  within  the  plate  from  the  gamma  photon  flux  incident  upon  it.     Since 
a  complete  energy  spectrum  of  gamma  photons  is  emitted  from  a  reactor,    it  is 
expedient,    for  purposes  of  calculation,    to  divide  the  spectrum  into  several 
discrete  energy  groups.     A  set  of  flux  data  taken  from  the  reactor  of  a  pro- 
totype nuclear  power  plant  was  employed  (1?).     The  spectrum  is  divided  into 
eight  energy  groups,    and  all  subsequent  calculations  are  based  on  these  data. 

For  a  steel  plate  at  a  point  1.5  meters  from  the  reactor,  the  data  pre- 
sented are  shown  in  Table  1 .  The  absorption  coefficients  for  both  steel  and 
aluminum  are  also  given  in  this  table. 

Table  1 .     Flux  data  from  a  prototype  nuclear  power  plant  and 

energy  absorption  coefficients  of  steel  and  aluminum 
at  these  correst»nding  energies. 


Ener 


Energy 
(Mev) 


Heat  flux 
(1000  Btu/hr,ft2) 


Energy  absorption  coefficient  (ft"1)* 


Steel 


Aluminum 


6.50 

3.05293 

4.72 

4.95276 

3.^0 

16.66158 

2.06 

24.02722 

1.26 

16.25520 

0.71 

13.51214 

0.41 

10.87067 

0.15 

5.28294 

5.5839 
5.4913 
5.3675 
5.5260 

5.9802 

6.7056 

7.0744 

19.4066 


1.5377 
1.5978 
1.7044 

1 .8983 
2.1269 
2.3180 

2. 361 3 
2.3287 


*  see  reference  (9). 

Since  the  number  of  mean  free  paths  through  any  of  the  steel  plates  con- 
sidered  is  small,  the  effect  of  buildup  will  be  minimal  (3).  Hence,  the  a- 


mount  of  Compton  scattering  is  small  and  the  linear  energy  absorption 


]k 


coefficient  gives  a  reasonably  accurate  description  of  the  attenuation  of 
the  gamma  photons  in  the  plate.   Therefore,  in  the  solution  to  Bq.(5)t  the 
first  step  is  to  sum  the  flux  and  absorption  coefficient  over  each  of  the 
eight  energy  groups.  It  is  then  possible  to  iterate  for  the  value  of  the 
equivalent  gray  body  connective  coefficient,  hr,  by  successively  solving 
1*1.(6),  (5)  and  (?).  This  procedure  is  outlined  in  greater  detail  in  Appen- 
dix D. 

A  complete  investigation  of  the  temperature  of  a  1/4,  1/2,  3/^t  and  1 
inch  plate,  from  zero  time  until  the  steady  state  temperature  distribution 
had  been  reached,  was  undertaken.   For  each  calculation  made  in  this  study 
the  base  temperature  was  assumed  to  be  500°R.   The  free  plate  thermal  stress 
distribution  was  also  calculated  at  various  times.  The  time  increment,  chosen 
was  0.01  hr.,  and  the  stipulation  was  made  that  this  increment  would  be 
doubled  each  time  the  inside  temperature  was  less  than  1 5°F  higher  than  that 
of  the  previous  time  considered  at  this  same  position.   Steady  state  was  con- 
sidered to  have  been  reached  when  the  inside  temperature  was  within  10°F  of 
the  steady  state  component  of  the  temperature  in  Bi.(5). 

Besserer  (2)  has  presented  data  on  the  physical  properties  of  steel. 
Many  of  these  properties  were  found  to  be  temperature  dependent,  and  hence, 
it  was  necessary  to  alter  the  value  of  these  physical  constants  as  the  tem- 
perature of  the  plate  increased.  Rather  than  fit  the  curves  which  Besserer 
presents,  accurately,  a  linear  approximation  was  employed.  The  only  devia- 
tion from  this  procedure  ocourred  in  the  case  of  the  modulus  of  elasticity 
which  was  approximated  by  two  connecting  linear  curves.  The  properties  con- 
sidered to  be  temperature  dependent  were  the  coefficient  of  thermal  expansion, 
the  thermal  conductivity,  the  emissivity,  the  densi  +  y  and  the  modulus-  of  e. 
lasticity.   The  density  was  found  from  volumetric  considerations  (19),  and 
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is  given  by 

fM-jL- _.  (21) 

where  A-  is  the  density  at  some  base  temperature  T0.   The  heat  capacity  and 
Poisson's  ratio  were  considered  to  be  non  temperature  dependent.   In  each 
case  the  temperature  chosen  for  alteration  of  the  constants  was  Ihe  average 
temperature  of  the  plate  at  the  time  of  interest. 

For  the  case  of  aluminum,  a  point  was  chosen  at  which  the  level  of  the 
flux  was  1/30  of  the  level  at  1.5  meters.   The  energy  distribution  at  such 
a  location  was  not  available;  consequently  the  previous  data  were  employed, 
with  the  intensity  of  each  energy  group  reduced  by  a  factor  of  JO. 

Since  the  absorption  coefficient  of  aluminum  is  significantly  less  than 
that  of  iron,  the  buildup  is  negligible  and  the  linear  energy  absorption 
coefficient  was  again  used.   For  aluminum,  the  only  thicknesses  considered 
were  one  half  and  one  inch.   Because  of  the  low  absorption  coefficient  and 
because  of  the  low  emissivity,  the  temperature  rose  at  a  much  slower  rate 
than  was  found  for  steel,  but  for  a  corresponding  flux  level  the  steady  state 
temperature  was  higher  in  aluminum.  Thus  the  first  time  increment  considered 
for  this  reduced  flux  is  0.1  hr.   It  was  also  found  that  each  of  the  tempera- 
ture dependent  physical  quantities  could  be  approximated  by  the  same  type 
curves  as  were  used  for  steel. 


The  Single  Plate  with  Time  Dependent  Gamma  Photon  Flux 
and  Radiative  Heat  Loss 


In  considering  the  temperature  distribution  in  a  metal  plate  when  it 
is  subjected  to  a  time  dependent  gamma  photon  flux,  it  is  convenient  to  em- 
ploy the  use  of  finite  differences.  Hence,  it  is  essential  that  the  size  of 
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the  time  increment  be  carefully  considered  before  attempting  to  find  a  solu- 
tion.  It  was  found  that  in  the  case  of  steel,  a  time  increment  of  0.01  hr. 
was  small  enough  to  give  answers  which  were  in  excellent  agreement  with  those 
obtained  from  the  analytic  calculations  for  a  1  jk   inch  plate.  Small  values 
of  the  ratio  At/ ax  tend  to  give  the  most  accurate  results  (18).   In  these 
calculations  ten  subdivisions  were  considered  for  all  plate  thicknesses. 
Hence,  the  aocuracy  of  the  calculations  for  a  1  /^  inch  plate  were  sufficient 
to  insure  good  results  for  thicker  plates. 

Maienschein,  et  al.  (13)  present  data  on  the  emission  of  gamma  photons 
and  their  energy  distribution  from  fission  products  formed,  after  various 
times  of  reactor  operation.  The  graphs  presented  show  that  the  amount  of 
higher  energy  radiation  present  decreased  very  rapidly  for  all  periods  of 
reactor  operation.  Since  there  is  a  higher  concentration  of  photons  of  a 
lower  energy  present,  it  would  be  possible  to  braoket  the  photons  of  all  en- 
ergies into  a  single  energy  group  at  the  average  energy,  provided  this  aver- 
age energy  is  olose  to  the  energy  peak  (?). 

The  average  energy  for  each  of  several  periods  of  reactor  operation  was 
calculated.  This  was  done  by  using  Simpson's  rule  (20)  for  energy  increments 
of  0.125  Mev.  in  the  range  0.25  to  5.50  Mev.  The  average  energy  for  the  fis- 
sion products  of  different  periods  of  reaotor  operation  varied  from  1.02  to 
1.29  Mev.,  with  the  fission  produots  of  longer  operating  time  having  the  low- 
er average  energy.  This  was  also  found  to  be  near  the  energy  peak,  especially 
for  the  longer  periods  of  reaotor  operation.  Hence,  the  use  of  a  single  en- 
ergy group  at  the  average  energy  for  the  purpose  of  calculation  is  justified. 

The  longest  operating  time  considered  in  the  above  data  is  about  25  min- 
utes. Operating  times  of  longer  than  about  12  minutes  yield  essentially  the 
same  average  energy  photons.  Sinoe  the  operating  times  considered  were  longer 
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than  12  minutes,    a  single  energy  of  1.05  Mev.  was  used  for  obtaining   the  ab- 
sorption coefficient.     Using  this  procedure,    it  was  also  more  expedient   to 
consider  the  total  energy  flux  in  preference  to  taking  a  summation  over  each 
of  the  eight  energy  groups. 

A  continuous  description  of  the  temperature  in  the  plate  was  obtained 
by  employing  the  analytic  expression  described  above,   until   the  time  of  re- 
aotor  shutdown.      The  period  of  reactor  operation  was  varied,    and  for  each  of 
the  several  oases  considered,    the  equations  in  finite  difference  form  were 
used  to  find  the  temperature  distribution  at  the  succeeding  time  increments. 
The  choice  of  the  point  of  termination  of  these  calculations  is  completely 
arbitrary  and  depends  on  how  closely  the  system  is  desired   to  be  investigated, 
at  each   succeeding  time  increment  the  flux  was  found  from  the  equations  pre- 
sented by  Lottes  (12)  and  inserted  into  the  matrix. 

As  has  been  previously  dJ3cussed,    the  temperature  dependent  constants 
in  the  matrix  must  be  ohangod   *fter  o»*h  iteration.     Since  no  analytic  expres- 
sion for  the  temperature  distribution  was  available,   it  was  not  possible  to 
directly  calculate  the  average  temperature.      Instead,    this  was  done  by  using 
Simpson's  rule  with  the  temperature  distribution  found  from  the  numerical 
solution  to  the  matrix. 


Two  Parallel  Plates  with  both  Convection  and  Radiation 
Heat  Loss  in  a  Constant  Gamma  Photon  Beam 


In  attempting  to  provide  additional  cooling  by  the  use  of  two  parallel 
plates  with  gas  flow  between  them,   only  the  steady  state  case  was  considered. 
It  was  felt  that  the  operation  of  the  reactor  for  lesser  periods  of  time  would 
not  raise  the  temperature  or  the  thermal  stress  levels  to  abnormal  proportions, 
The  investigation  of  the  effect  of  a  number  of  convection  coefficients  on  the 
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magnitude  of  the  temperature  in  steel  plates  and  the  distribution  of  the  two 
plate  thicknesses  necessary  to  keep  the  average  temperature  of  the  plates 
equal,  was  made.   The  range  of  convection  coefficients  used  was  from  1  to 
25  Btu/hr,ft2,OF. 

It  was  again  necessary  to  alter  the  temperature  dependent  physical  prop- 
erties of  each  of  the  two  plates  as  the  average  temperatures  changed.  How- 
ever, it  was  noted  that  the  sign  of  some  of  these  properties  changed  if  the 
temperature  became  too  high.   This  phenomenon  occurs  because  the  linear  ap- 
proximations used  are  invalid  at  excessively  high  temperatures,  especially 
near  the  melting  point.   Since  some  original  division  of  the  plate  thickness 
was  necessary,  it  was  possible  that  for  low  convection  coefficients  the  tem- 
perature of  the  inside  plate  would  be  abnormally  high  until  its  thickness  had 
been  reduced  enough  to  lower  the  temperature  to  a  reasonable  level.   There- 
fore, some  of  the  physical  constants  would  not  be  represented  by  their  true 
values.   In  general  this  would  not  hinder  the  iteration.  However,  as  the 
thickness  of  the  inside  plate  is  reduced,  there  exists  a  point  at  which  a 
change  in  the  sign  of  the  physical  constants  occurs.  In  the  case  of  the  ther. 
mal  conductivity,  one  segment  of  the  equation  involved  in  the  calculation 
of  the  temperature,  contains  its  reciprocal  so  that  when  it  arrives  at  the 
point  where  the  sign  ohanges,  the  value  of  this  portion  of  the  equation  var- 
ies from  a  high  positive  to  a  high  negative  number.   Consequently,  this  iter- 
ation would  probably  never  converge.   If  the  average  temperature  of  the  plate 
reached  this  unrealistic  level,  it  was  set  equal  to  the  steady  state  average 
temperature  of  a  single  plate  of  equivalent  thickness,  for  the  purpose  of 
altering  the  physical  constant  involved.   From  this  point  the  calculations 
proceeded  in  a  nonnal  manner  until  the  average  temperatures  of  the  two  plates 
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were  within  the  Units  specified.     The  limits  set  were  that  the  temperature 
difference  between  the  plates,  divided  by  their  combined  average  temperature, 
be  less  than  1<H*. 
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DISCUSSION  AND  RESULTS 


The  Single  Plate  with  Constant  Gamma  Photon  Flux  and 
Radiative  Heat  Loss 


Figures  1  and  2  show  the  theoretically  predicted  temperature  responses 
of  two  different  sized  plates  of  steel  and  aluminum  to  a  constant  gamma  pho- 
ton flux.  The  calculations  on  the  aluminum  plate  were  made  with  the  flux 
level  of  each  of  the  eight  energy  groups  being  1  /30  of  the  flux  level  to  which 
the  steel  plate  was  exposed. 

The  temperature  rise  in  both  the  1  fk   inch  and  the  one  inch  steel  plates 
shown  in  Fig.  1  is  quite  rapid.  In  the  1/4  inch  plate  the  average  heat  gen- 
eration rate  is  somewhat  higher  and  the  initial  temperature  rise  is  greater 
than  in  the  one  inch  plate.  However,  after  a  short  period  of  time,  the  thin- 
ner plate  is  able  to  conduct  a  greater  percentage  of  the  heat  generated  with- 
in it  to  the  surface.   Consequently,  its  rate  of  ascension  quickly  begins  to 
fall  while  the  temperature  in  the  one  inch  plate  continues  to  rise  in  almost 
a  linear  manner.   Phis  is  also  true  because  of  the  lesser  amount  of  total 
heat  capable  of  being  stored  within  a  plate  of  smaller  dimensions. 

An  indication  of  the  effect  of  the  thermal  conductivity  on  the  rate  of 
heat  removal  can  be  obtained  by  examining  the  change  in  the  temperature  dis- 
tribution in  each  of  the  two  plates  as  the  period  of  exposure  to  the  gamma 
flux  progresses.  This  is  done  in  Tables  2  and  3.  The  time  of  the  first  dis- 
tribution shown  in  each  case,  occurs  during  the  period  when  there  is  a  linear 
temperature  rise  and  the  last  distribution  is  at  steady  state.  At  first, 
both  the  plates  move  slowly  from  the  initially  flat  temperature  distribution. 
However,  the  temperature  differential  between  the  inside  and  outside  tempera- 
tures in  the  one  inch  plate  increases  quite  rapidly  compared  to  the 
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Table  2.  Temperature  distribution  in  a  1/4  inch  steel  plate  after 

irradiation  by  a  constant  gamma  photon  beam  for  various  times 


Distance  through 
the  plate     (ft) 

Temperature 

Time  of  reactor 

no 

operation     (hr) 

0.020833 

601.818 

0.01 

0.018750 

601 .827 

0.01666? 

601  .842 

0.014583 

601.859 

0.012500 

601  .880 

0.010417 

601  .890 

0.008333 

601  .922 

0.006250 

601  .941 

0.004167 

601 .955 

0.002083 

601 .966 

0 

601 .969 

0.020833 

1344.598 

0.10 

0.01 8750 

1344.972 

0.01 6667 

1345.315 

0.014583 

1345.624 

0.012500 

1345.897 

0.010417 

1346.134 

0.008333 

1346.331 

0.006250 

1346.489 

0.004167 

1346.603 

0.002083 

1346.673 

0 

1346.697 

0.020833 

1707.216 

0.34 

0.01 8750 

1708.341 

0.01 6667 

1709.357 

0.014583 

1710.261 

0.012500 

1711.051 

0.010417 

1711.725 

0.008333 

1712.282 

0.006250 

1712.719 

0.004167 

1713.033 

0.002083 

1713.224 

0 

1713.288 

0.020833 

1734.257 

steady 

0.01 8750 

1735.469 

state 

0.016667 

1736.563 

0.014583 

1737.536 

0.012500 

1738.385 

0.010407 

1739.110 

0.008333 

1739. 708 

0.006250 

1740.177 

0.004167 

1740.515 

0.002083 

1740.719 

0 

1740.788 
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Table  3«   Temperature  distr.ibu  tior  in  a  one  inch  steel  plate  after 

irradiation  by  a  cor?* art  gamma  photon  beam  for  various  times 


Distance  through 
the  plate  (ft) 


Temperature 
(°R) 


Time  of  reactor 
operation     (hr) 


0.083333 
0.075000 
O.O66667 
O.058333 
0.050000 
0.041667 

0.033333 
0.025000 
0.01 6667 
0.008333 
0 

O.O83333 
0.075000 

O.066667 

0.058333 
0.050000 

0.041667 

0.033333 
0.025000 

0.01 6667 

0.008333 

0 

0.083333 
0.075000 
0.066667 
0.058333 
0.050000 
0.041667 
0.033333 
0.025000 
0.01 6667 

0.008333 
0 

0.083333 
0.075000 
0.066667 

0.058333 
0.050000 

0.041667 

0.033333 
0.025000 

0.01  6667 

0.008333 
0 


579.270 
579.429 
579.834 
580  >36 
581. 178 
582.006 

582.853 
583.649 
584.327 
584.797 
584.974 

1  291  .  625 
1293.035 
1294.622 

1296.325 
1298.076 
1299.804 
1301  .430 
1302.868 

1304.027 
1304.805 
1305.090 

21  31  .  287 
2145.108 
2157.941 
21 69. 698 
2180.286 
2189.598 
2197.524 
2203.941 
2208.714 
2211 .696 
2212.729 

2319.233 
2340.928 
2360.878 
2378.980 
2395.123 
2409.188 
2421  .043 
2430.550 

2437.555 
2441 .893 
2443.382 


0.01 


0.10 


0.34 


~7^ 


steady 
state 
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differential  for  the  1  /k   inch  plate.  A  more  detailed  description  of  this 
response  is  shown  in  Fig.  3. 

For  the  case  of  the  aluminum  plates  shown  in  Fig.  2,  the  temperature 
rise  is  considerably  slower  than  that  of  steel.  In  fact,  for  the  magnitude 
of  the  flux  considered,  the  time  necessary  to  approach  steady  state  is  about 
ten  times  as  great  for  the  aluminum  plate  as  for  the  steel  plate.  Calcula- 
tions show  that  when  the  aluminum  plate  is  subjected  to  the  same  gamma  flux 
as  the  steel  plate,  it  reaches  temperatures  of  several  hundred  degrees  Fahr- 
enheit higher  than  the  steel  plate  for  all  the  sizes  of  plates  examined. 
However,  the  temperature  rise  remains  considerably  slower  in  the  aluminum. 

This  phenomenon  is  a  result  of  several  of  the  physical  properties  in- 
volved. Because  the  aluminum  has  a  significantly  lower  absorption  cross 
section  than  steel  (see  Table  1),  the  amount  of  heat  generated  within  a  com- 
parable volume  of  aluminum  is  distinctly  less  than  that  generated  in  steel 
during  any  given  time  increment.   Consequently,  the  temperature  has  a  tend- 
ency to  ascend  at  a  much  slower  rate  in  aluminum.  However,  the  temperature 
rise  is  also  affected  by  the  thermal  conductivity  and  the  emissivlty  of  the 
material  involved.  The  high  thermal  conductivity  of  aluminum  allows  the  heat 
generated  within  the  plate  to  reach  the  surface  far  more  rapidly  than  is  the 
case  for  steel.  But  since  the  emissivity  of  aluminum  is  low,  the  heat  has 
a  reduced  chance  to  be  radiated  once  it  reaches  the  surface. 

By  far,  the  overriding  difference  between  aluminum  and  steel  during  the 
early  portion  of  the  heating  period  is  the  low  absorption  cross  section  of 
aluminum.   The  amount  of  heat  generated  is  not  nearly  enough  to  cause  as  rap- 
id a  temperature  rise  as  is  found  in  steel.  However,  because  of  the  low  em- 
issivity of  aluminum,  it  retains  a  significantly  greater  percentage  of  the 
heat  generated  within  it.   Consequently,  when  the  aluminum  plate  finally 
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reaches  an  equilibrium  condition,  the  total  amount  of  heat  retained  within 
it  is  greater  than  for  similar  conditions  in  steel.  Thus,  the  temperature 
in  aluminum  reaches  a  higher  steady  state  value. 

It  is  also  seen  from  Fig.  2,  that  for  all  practical  purposes,  the  temp- 
erature of  both  the  one  half  inch  and  one  inch  aluminum  plates  is  the  same 
for  a  period  of  two  hours.   This  can  be  attributed  in  part  to  the  high  ther- 
mal conductivity  of  the  aluminum,  which  allows  a  more  rapid  flow  of  heat  a- 
cross  any  thickness  of  plate.   The  significance  of  this  high  conductivity 
is  attested  by  the  fact  that  the  temperature  drop  across  a  one  inch  aluminum 
plate  at  the  steady,  state  temperature  distribution  is  less  than  one  degree 
Fahrenheit,  while  for  steel  this  temperature  differential  exceeds  100°F. 

This  high  conductivity  of  the  aluminum  approximates  one  of  the  stipula- 
tions for  the  "thin  body"  temperature  transients  discussed  by  Unterberg  (22). 
Such  a  system  assumes  zero  internal  resistance  to  heat  flow,  so  that  no  temp- 
erature gradient  exists  within  the  material.   This  assumption  is  satisfied 
by  an  infinite  (k/x),  and  removes  the  spatial  temperature  variation  from  the 
system.  Hence,  the  heat  balance  for  the  case  of  radiation  heat  loss  is  given 

by 

V^^X=HV-Ao-e(T^x<)  (22) 

where  V  is  the  volume  of  the  body  and  A  is  the  area  of  the  surface  radiating 
to  the  surroundings. 

The  steady  state  case  is  then  given  by  (dT/dt)  =  0,  whence 
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where 


(2<0 


7  =  fV-e.  (25) 

If  a  quantity,  p,  is  defined  as 


1-fiB 


(26) 


the  temperature  is  given  by  the  transcendental  equation 

2fUt  =  ^CfT)  -^CfT)  (27) 

where 

^(fT)   .1wK-'(-pT)  +  W(fT) 


(28) 


Solution  of  Bq.(23)  for  a  one  half  inch  aluminum  plate,  shows  the  steady 
state  "thin  body"  temperature  to  be  1 337°R.   For  a  one  inch  plate  this  temp- 
erature is  1  553°R.   The  exact  solution  for  these  temperatures,  shows  them 
to  be  1026°R  and  1 1  ?3°R  respectively.   The  temperature  rise  for  the  Mthin 
body"  aluminum  takes  on  the  same  relative  Kiiape,  but  as  the  temperature  in- 
creases, the  departure  from  the  exact  solution  becomes  more  apparent.   Thus, 
the  "thin  body"  approximation  for  the  aluminum  plates  examined  starts  to  ap- 
proach the  exact  solution,  but  does  not  give  accurate  answers  except  at  temp- 
eratures only  slightly  above  the  base  temperature. 

It  is  possible  to  consider  the  thermal  stresses  expressed  in  Eq.(8)  by 
using  the  equation  in  its  entirety  or  by  considering  any  of  its  several  parts, 
The  simplest  case  is  to  examine  the  first  term  of  Eq.(8)  alone.   This  is 
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equivalent  to  holding  the  ends  of  the  plate  in  perfect  restraint,  and  gives 
the  maximum  thermal  stress  to  which  the  plate  may  be  subjected  (21).   Al- 
though this  condition  could  probably  not  be  precisely  reconstructed  in  a 
nuclear  rocket,  a  sketch  of  the  physical  representation  is  given  in  Fig.  3a. 
Here,  the  plate  is  pinned  at  the  two  ends  to  supports  which  remain  fixed. 
Consequently,  the  plate  can  neither  expand  longitudinally  nor  bend  about 
either  end  to  relieve  the  applied  stresses.   If  the  physical  constants  in- 
volved in  this  calculation  are  assumed  to  be  constant,  the  maximum  thermal 
stress  is  a  linear  function  of  the  temperature. 

Figure  ^  shows  the  maximum  yield  stress  for  ^1  30  steel.   The  calculation 
of  the  maximum  thermal  stress  for  steel  shows  it  to  exceed  the  maximum  yield 
stress  by  the  time  the  temperature  has  risen  more  than  200°F  above  the  base 
temperature.   For  the  one  inch  plate,  which  has  the  slowest  initial  tempera- 
ture rise  of  any  of  the  cases  investigated,  the  time  necessary  to  attain  this 
temperature  is  less  than  two  minutes.  Since  it  is  highly  illogical  that  any 
nuclear  rocket  power  plant  would  feasibly  remain  in  operation  for  that  short 
a  period  of  time,  provisions  must  obviously  be  made  to  allow  for  expansion 
of  the  system  involved.   Since  the  maximum  thermal  stress  attains  a  value 
of  five  times  the  yield  stress  within  ten  minutes,  this  consideration  is  ex- 
tremely important. 

In  examining  the  maximum  thermal  stress  in  aluminum,  the  results  are 
found  to  be  considerably  different  than  for  steel,  as  is  shown  in  Fig.  5» 
The  point  at  which  the  maximum  thermal  stress  first  exoeeds  the  yield  stress 
for  7075-T6  aluminum  with  the  flux  level  considered,  does  not  occur  until 
almost  two  hours  have  elapsed.  Thus  for  periods  of  reactor  operation  in  the 
1  5  to  30  minute  range,  the  thermal  stresses  in  aluminum  when  subjected  to 
the  heating  level  discussed,  are  well  within  the  destructive  limits. 
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FIG.    3a!     SKETCH    OF    FLAT     PLATE    HELD    IN     PERFECT 
RESTRAINT. 
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FIG.    3b!      SKETCH    OF    FLAT     PLATE     RESTRAINED 
ONLY       IN     BENDING. 
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FIG.    3c:     CANTILEVER    ARRANGEMENT     IN     WHICH 

PLATE     MAY     MOVE     IN     ANY     DIRECTION. 
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It  is  also  possible  to  consider  the  thermal  stress  produced  in  a  plate 
which  is  restrained  in  bending.  This  is  done  by  evaluating  the  first  two 
terms  in  Eq.(8).  Further  examination  shows  that  this  equation  is  simply 

where  Q   represents  the  average  temperature  differential  above  the  base 

SIV 

temperature.  This  case  is  represented  schematically  in  Fig.  3b.  This  is  a 
cantilever  type  arrangement  in  which  the  unattached  end  of  the  plate  is  free 
to  expand  longitudinally,  but  is  held  between  two  immovable  guides  so  that 
it  may  not  bend  about  its  fixed  end.  The  quantity  which  is  of  primary  inter- 
est here  ia  again  the  maximum  possible  thermal  stress.  This  is  found  by  con- 
sidering the  point  at  which  the  largest  departure  from  the  average  tempera- 
ture exists. 

Figure  6  shows  the  maximum  thermal  stress  of  several  different  sizes 
of  steel  plates  which  are  restrained  in  bending.   For  each  of  these  plates 
the  maximum  thermal  stress  exists  at  the  outer  extremity  of  the  plate.   Also, 
the  stress  is  always  found  to  be  one  of  tension  at  this  point  in  the  plate. 
It  is  seen  that  as  the  plate  becomes  thicker,  the  maximum  thermal  stress  in- 
creases at  a  rapid  rate,  especially  at  higher  temperatures.  Thus,  in  the 
case  of  thermal  stresses  there  is  a  definite  size  of  plate  which  may  be  used 
if  it  is  desired  to  keep  the  thermal  stresses  in  the  range  which  is  well  be. 
low  the  yield  stress  of  steel. 

It  should  be  noted  in  Fig.  4,  however,  that  the  maximum  thermal  stress 
for  a  plate  restrained  in  bending  does  not  begin  to  seriously  approach  the 
yield  stress  in  steel  until  the  temperature  of  the  plate  has  risen  above 
2000°R.   For  the  flux  levels  which  were  considered  here,  this  case  does  not 
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arise  until  the  thickness  of  the  plate  is  approximately  one  half  inch.   For 
plates  of  this  and  greater  thicknesses,  a  temperature  of  2000°R  is  reached 
in  approximately  15  minutes.   For  plates  of  the  thicknesses  examined,  this 
maximum  thermal  stress  never  exceeds  the  yield  stress.  wowever,  plates  of 
larger  thickness  would  undoubtedly  tend  to  have  their  maximum  thermal  stresses 
increase  at  an  increasing  rate.  This  is  evidenced  by  the  rate  of  increase 
of  the  temperature  differential  across  the  plate  as  the  thickness  is  increas- 
ed (see  Fig.  3).  This  rate  of  increase  can  also  be  seen  by  examining  the 
steady  state  line  in  Fig.  6,  which  gives  a  more  detailed  description  of  the 
maximum  thermal  stress  for  a  plate  restrained  in  bending  for  each  of  the 
thicknesses  considered. 

In  Fig.  k,    the  steady  state  line  has  been  extrapolated  slightly  beyond 
the  calculated  data  in  an  attempt  to  give  some  indication  of  the  maximum 
thermal  stress  of  a  plate  larger  than  those  discussed  here.  However,  it 
must  be  remembered  at  this  point,  that  one  of  the  assumptions  used  in  the 
calculations  of  this  study  was  that  the  plates  under  consideration  were  thin 
enough  so  that  Compton  scattering  could  be  neglected.  This  assumption  be- 
gins to  rapidly  break  down,  especially  for  the  lower  energy  groups,  if  cal- 
culations are  made  on  thicker  plates.  When  the  amount  of  Compton  scattering 
begins  to  increase  significantly,  the  total  energy  release  predicted  by  the 
equations  derived  in  this  study  is  seriously  underestimated.  Applying  this 
information  to  the  thermal  stresses  under  consideration,  the  magnitude  of 
these  thermal  stresses  would  increase  at  an  even  more  rapid  rate  than  is 
predicted  by  the  equations  presently  employed. 

In  the  case  of  aluminum,  the  temperature  differential  across  the  plate 
was  found  to  be  almost  non  existent.  As  previously  discussed,  this  is  due 
to  the  high  thermal  conduct! /Lty  of  aluminum.   Consequently,  the  maximum 
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thermal  stress  for  either  of  the  plates  considered,  when  restrained  in  bend- 
ing, lies  in  the  range  of  zero  to  ten  psi.   This  can  be  considered  as  neg- 
ligible for  all  practical  purposes.   The  ratio  of  the  maximum  thermal  stress 
for  an  aluminum  plate  restrained  in  bending  to  the  yield  stress  is  1 0~J   or 
less,  so  that  the  plate  essentially  does  not  "feel"  the  thermal  stress  which 
exists  within  it.   The  "thin  body"  approximation  then  seems  to  gain  a  slight 
degree  of  stature  since  it  predicts  a  zero  stress  in  all  cases  except  that 
of  perfect  restraint. 

The  final  manner  in  which  the  thermal  stress  may  be  examined  is  to  con- 
sider a  free  plate.   In  this  case  Sq.(8)  is  used  in  its  entirety.   The  rep- 
resentation of  such  a  physical  construction  is  shown  in  Fig.  3C«  The  one 
end  of  the  plate  is  fixed  with  the  plate  free  to  expand  in  any  direction 
from  this  end.   Such  a  system  is  evidenced  in  any  simple  cantilever  arrange- 
ment.  Figure  7  gives  an  indication  of  the  thermal  stress  response  of  a  1 \h 
inch  steel  plate.   In  attempting  to  explain  this  unusual  transient  response, 
it  is  helpful  to  examine  the  individual  terms  which  determine  the  free  plate 
stress. 

The  second  term  in  this  equation  is  simply  a  constant.   The  last  term 
is  linear  with  position,  and  is  zero  at  the  midpoint  of  the  plate.   The  first 
term  is  simply  a  constant  times  the  negative  of  the  temperature  distribution. 
Figure  8  shows  the  temperature  distribution  at  three  different  times.   After 
the  higher  order  transient  terms  of  Eq.(5)  have  died  out,  the  temperature 
profile  begins  to  closely  approximate  a  cosine  distribution.  The  free  plate 
thermal  stress  for  such  a  distribution  has  the  appearance  of  the  curve  shown 
for  t  =  0.12  hr.  in  Fig.  7. 

Examination  of  the  temperature  profile  for  t  =  0.01  hr.  shows  an  in- 
flection point.   Consequently,  this  temperature  profile  has  the  appearance 
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of  nearly  180°  of  a  cosine  curve.   The  resultant  thermal  stress  also  extends 
through  two  quadrants  of  the  normal  cosine- type  distribution.  Within  a  very 
short  time  inflection  point  disappears  from  the  temperature  profile  and  the 
free  plate  thermal  stress  distribution  begins  to  again  take  on  the  appear- 
ance of  a  simple  cosine- type  distribution.   Although  these  curves  are  rather 
interesting  to  investigate,  they  do  not  present  any  serious  physical  diffi- 
culties with  regard  to  structural  considerations. 

The  one  inch  plate  exhibits  the  same  cosine  type  distribution  as  the 
1  /4  inch  plate,  but  it  takes  slightly  longer  before  the  higher  order  temp- 
erature transient  terms  have  been  sufficiently  damped  to  cause  the  thermal 
stress  curve  to  show  a  simple  cosine  determined  distribution.   The  magnitude 
of  the  maximum  thermal  stress  in  a  free  plate  is  about  2000  psi  for  a  one 
inch  steel  plate  at  steady  state.   In  no  case  does  the  maximum  free  plate 
stress  exceed  that  of  a  plate  which  is  completely  restrained  in  bending. 
Thus,  this  former  stress  should  be  of  little  consequence  in  considering  the 
structural  feasibility  of  any  rocket  construction. 

As  would  be  expected  from  the  case  of  the  steel  plate  just  discussed, 
the  maximum  free  plate  thermal  stress  in  aluminum  is  also  less  than  the 
maximum  thermal  stress  for  a  plate  restrained  in  bending.   The  free  plate 
stresses  exhibit  the  same  type  distribution  as  was  evidenced  in  steel,  how- 
ever, these  are  so  small  that  they  can  be  completely  neglected.  The  maxi- 
mum thermal  stress  for  any  of  the  cases  examined  in  aluminum  was  found  to 
be  less  than  one  psi. 


The  Single  Plate  with  Time  Dependent  Gamma  Photon  Flux 
and  Radiative  Heat  Loss 


It  is  to  be  expected  that  a  time  dependent  investigation  of  the 


no 


0.1      0.2 


0.24 


025 


0.21  0.22  0.23 

TIME,     t,    (HR) 

FIG.  9.        OUTSIDE     TEMPERATURE    OF    1/2    INOH    STEEL 
PLATE     AFTER      SHUTDOWN     WHEN    SUBJECTED  TO   TIME 
DEPENDENT    GAMMA      PHOTON    FLUX 


*n 


temperature  within  a  flat  plate  will  present  a  smooth  curve.  Figure  9  shows 
the  curve  for  the  case  of  shutdown  after  12  minutes  of  reactor  operation. 
The  time  scale  after  shutdown  has  been  expanded  in  an  attempt  to  give  a  more 
distinct  appraisal  of  the  temperature  during  this  period.  Because  of  the 
large  temperature  increments  necessary  to  plot  this  graph,  it  is  difficult 
to  distinguish  that  for  a  short  period  the  curve  is  concave  downward,  after 
which  it  goes  through  an  inflection  point  and  reverses  itself.   For  the  time 
increments  chosen,  the  curve  appears  to  show  a  discontinuity.  Table  ^  gives 
che  actual  temperature  distribution  immediately  before  and  0.001  hr.  after 
shutdown  for  two  different  times  of  reactor  operation.   It  is  most  probable 
that  this  apparent  discontinuity  would  be  eliminated  by  choosing  a  time  in- 
crement small  enough.   However,  the  increment  could  not  be  shortened  below 
0.1  seconds  since  this  is  the  endpoint  of  the  graph  for  the  time  dependent 
flux  equation  given  by  Lottes  (12). 

A  comparison  was  also  made  of  the  temperature  of  a  steel  plate  using 
the  time  dependent  flux  equation  of  Lottes  and  also  assuming  zero  flux  im- 
mediately after  shutdown.   Tables  5  and  6  show  these  temperatures  for  several 
different  times  after  shutdown.   The  operating  time  of  the  reactor  was  taken 
as  about  22  minutes,  at  which  time  the  rate  of  increase  of  the  temperature 
has  slowed  considerably  from  the  original  rate  of  ascension.   Examination 
of  these  tables  shows  that  after  a  period  of  two  to  three  minutes  the  dif- 
ference in  the  temperatures  between  the  two  cases  is  about  five  degrees  Fahr- 
enheit. 

Thus  the  amount  of  heating  provided  by  the  gamma  photons  after  shutdown 
is  quite  minimal  when  compared  to  the  heating  during  reactor  operation.  For 
the  case  of  zero  flux  after  shutdown,  the  temperature  would  asymptotically 
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Table  4.      Temperature 

distribution  in  a  1/2  in< 

3h   steel   plate  at  shut- 

down  and   0.001    hours  after   shutdown  for  operating   times  of 

0.20  and   0.37  hours 

Distance   through 

Temperature 

Time 

the  plate      (ft) 

(°R) 

(hr) 

0.0^1667                                       1779.283                                       0.200 

0.037500                                       1782.052 

0.033333                                     1784.592 

0.0291 67                                     1786.888 

0.025000                                     1788.929 

0.020833                                     1790.699 

0.01 6667                                     1792.184 

0.012500                                     1793.368 

0.008333                                     179^.236 

0.004167                                     I794.77O 

0                               179^.952 

0.041667                                   1773.911                                   0.200 

0.037500                                   1776.7^9                                      + 

0.033333                               1779.271                                0.001 

0.029167                                   1751.490 

0.025000                                   1783.409 

0.020833                                  1785.030 

0.016667                                   1786.354 

0.012500                                   1787.379 

O.OO8333                                1788.104 

0.004167                                  1788.535 

0                                1788.675 

0.041667                                  1963.272                                  0.370 

0.037500                                  1967.730 

0.033333                                  1971.786 

0.029167                                  1975. ^26 

0.025000                                  1978.636 

0.020833                                  1981.400 

0.01 6667                         1983.702 

0.012500                                  1985.524 

0.008333                                  1986.849 

0.004167                                  1987.659 

0                                      1987.933 

0.041667                                195^.915                                0.370 

0.037500                                1959.^71                                    + 

0.033333                            1963.526                            0.001 

0.0291 67                                 1967.095 

0.025000                                 1970.1 85 

0.020833                                 1972.798 

0.016667                                 1 97J4.933 

0.012500                                 1976.590 

0.008333                                1977.767 

0.004167                                 1978.467 

0 

1978.698 
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Table  5»  Temperature  distribution  in  a  1/2  inch  steel  plate  for  various 
times  after  reactor  shutdown  at  0.37  hours  without  additional 
heating  after  shutdown 


Distance  through 

Temperature 

Time  after 

the  plate      (ft) 

(°R) 

shutdown  (hr) 

0.04166? 

1954.616 

0.001 

0.037500 

1959.168 

0.033333 

1963.219 

0.0291  67 

1966.783 

0.025000 

1969.867 

0.020833 

1972.474 

0.016667 

1974.604 

0.012500 

1976.255 

0.008333 

1977.428 

0.004167 

1978.125 

0 

1978.354 

0.041667 

1881  .357 

0.010 

0.037500 

1 885.1 30 

0.033333 

1 888 .  51  0 

0.0291  67 

1891  .502 

0.025000 

1894.104 

0.020833 

1896.312 

0.01 6667 

1898.1  23 

0.012500 

1899.533 

O.OO8333 

1900.539 

0.004167 

1901  .141 

0 

1901.340 

0 . 041 667 

1750.670 

0.030 

0.037500 

1753.335 

0.033333 

1755.728 

0.029167 

1757.843 

0.025000 

1 7 39. 677 

0.020833 

1761 .230 

0.01  6667 

1762.502 

0.012500 

1763.493 

0.008333 

1764.202 

0.004167 

1764.629 

0 

1764.773 

0.041667 

1649.840 

0.050 

0.037500 

1651.851 

0.033333 

1 653.657 

0.029167 

1655.251 

0.025000 

1 656.633 

0.020833 

1657.804 

0.016667 

1658.762 

0:012500 

1 659.508 

0.008333 

1660.042 

0.004167 

I66O.363 

0 

1660.472 
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Table  6.      Temperature  distribution  in  a  1/2  inch   steel  plate  for  various 
times  after  reactor   shutdown  at  0.37  hours  when  exposed   to  a 
time  dependent  gamma  photon  beam 


Distance   through 

Temperature 

the  plate      (ft) 

(°R) 

0 . 041  667 

1954.915 

0.037500 

1959.471 

0.033333 

1963.526 

0.029167 

1967.095 

0.025000 

1970. 185 

0.020833 

1972.798 

0.016667 

1974.933 

0.012500 

1976.590 

O.OO8333 

^977.767 

0.004167 

1978.467 

0 

1978.698 

0.041667 

1 883.638 

0.037500 

1 887.434 

0.033333 

1890.837 

0.029167 

1 893.850 

0.025000 

1 896.470 

0.020833 

1898.693 

0.01 6667 

1900.516 

0.012500 

1901.935 

0.008333 

1902.950 

0.004167 

1903.558 

0 

1903.759 

0.041667 

1754.832 

0.037500 

1757.528 

0.033333 

1759.947 

0.029167 

1762.085 

0.025000 

1763.941 

0.020833 

1765.514 

0.016667 

1766.802 

0.012500 

1767.806 

0.008333 

1768.524 

0.004167 

1768.957 

0 

1769.102 

0 . 041  667 

1655.241 

0.037500 

1 657.283 

0.033333 

1659.111 

0.029167 

1 660.729 

0.025000 

1662.135 

0.020833 

1663.326 

0.01 6667 

1664.302 

0.012500 

1665.062 

0.008333 

1665.605 

0.004167 

1665.930 

0 

1666.037 

Time  after 
shutdown   (hr) 


0.001 


0.010 


0.030 


0.050 
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approach  the  base  temperature  of  the  plate.   If  a  plate  were  subjected  to  a 
time  dependent  gamma  flux  such  as  occurred  after  shutdown  of  a  reactor,  it 
would  show  some  transient  temperature  response.   For  the  case  examined,  where 
the  plate  is  first  heated  by  a  constant  gamma  flux  during  reactor  operation, 
the  transient  temperature  analysis  will  first  show  a  rapid  cooling  period. 
However,  the  temperature  will  eventually  approach  the  values  predicted  for 
a  plate  which  had  not  been  subjected  to  the  gamma  flux  during  reactor  opera- 
tion.  After  a  period  of  one  or  two  hours  this  flux  would  not  be  significant 
enough  to  keep  the  temperature  of  the  plate  much  above  the  base  temperature. 


Two  Parallel  Plates  with  both  Convection  and  Radiation 
Heat  Loss  in  a  Constant  Gamma  Photon  Beam 


In  attempting  to  minimize  the  thermal  stresses  in  each  of  two  parallel 
plates  which  have  a  gas  flowing  between  them,  the  minimization  of  the  max- 
imum thermal  stress  of  a  plate  perfectly  restrained  at  the  ends  was  chosen. 
This  is  accomplished  by  simply  equating  the  average  temperature  in  each  of 
the  two  plates.  However,  it  does  not  necessarily  follow  that  the  maximum 
thermal  stress  of  a  free  plate  or  of  a  plate  restrained  only  in  bending  will 
also  be  minimized. 

In  all  cases  examined  for  both  a  steel  and  an  aluminum  plate,  only  the 
steady  state  temperatures  and  thermal  stresses  were  calculated.   Rather  than 
seeking  the  transient  nature  of  the  double  plate  problem,  an  attempt  was 
made  to  correlate  the  resultant  temperature  and  thermal  stresses  to  the  mag- 
nitude of  the  convection  coefficients  used. 

Figure  1 0  shows  the  relationship  for  the  outside  steel  plate,  steady 
state  temperature  to  the  convection  coefficient.  When  the  convection  coef- 
ficient reaches  a  value  of  zero,  the  temperature  will  be  the  same  as  was 
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found  for  a  single  plate  of  the  same  total  thickness.  Because  this  condition 
would  require  the  reduction  of  one  of  the  plates  to  zero  thickness,  this 
case  could  not  be  directly  calculated.  However,  Fig.  10  has  been  extrapo- 
lated to  the  actual  single  plate  steady  state  temperature  and  the  curve  ap- 
pears to  be  smooth,  thus  giving  at  least  some  partial  verification  to  the 
calculations  performed.   The  calculations  made  for  aluminum  show  similar  type 
curves  which  smoothly  approach  the  single  plate  steady  state  temperature, 
as  is  shown  in  Fig .11. 

Figure  1 2  shows  the  percentage  of  the  total  steel  plate  thickness  which 
is  contained  in  the  outside  plate.   For  a  convection  coefficient  of  zero, 
each  plate  examined  will  have  its  entire  thickness  in  the  outside  plate. 
Consequently,  these  curves  were  also  extrapolated  to  that  value. 

Figure  13  shows  the  same  type  data  just  described,  when  obtained  for 
an  aluminum  plate.  The  flux  level  chosen  to  compute  this  data  was  1 /1 0  the 
flux  level  to  which  the  steel  plate  was  exposed.  This  level  is  quite  arbi- 
trary and  was  chosen  simply  to  keep  the  temperature  below  the  melting  point 
of  aluminum.  It  can  be  seen  that  the  relationship  of  the  plates  of  differ- 
ent thickness  to  each  other  on  this  graph  are  essentially  the  same  as  were 
found  for  steel.  However,  the  graphs  are  quite  dissimilar  in  other  respects. 

The  two  outstanding  differences  between  aluminum  and  steel  are  that  the 
division  of  the  aluminum  plate  is  much  more  nearly  equal  than  the  steel 
plate,  and  furthermore  retains  this  near  equality  until  the  convection  co- 
efficient becomes  almost  zero.   Both  these  differences  can  be  attributed  to 
the'  combined  effect  of  the  high  conductivity  of  aluminum  and  the  availabil- 
ity of  an  additional  mode  of  heat  transfer. 

In  considering  the  aluminum  plate  with  radiation  heat  loss  only,  the 
low  emissivity  of  the  aluminum  was  found  to  retain  the  heat,  despite  the 
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FIG.  II         EFFECT      OF     CONVECTION     COEFFICIENT     ON 
AVERAGE     STEADY     STATE     TEMPERATURE     IN    AN 

ALUMINUM      PLATE 
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CONVECTION     COEFFICIENT,    h,    (BTU/hr,  Ft.2»F) 

FIG.  12.      PERCENTAGE    OF  TOTAL  THICKNESS  IN  OUTSIDE   PLATE  AT  DIFFERENT 
CONVECTION    COEFFICIENTS    WHEN  AVERAGE  STEEL  PLATE   TEMPERATURES 

ARE      EQUAL 
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CONVECTION    COEFFICIENT,   h,    (BTU/hr.,  Ftf  *F) 
FIG.  13.      PERCENTAGE    OF     TOTAL    THICKNESS    IN    OUTSIDE   PLATE 
AT  DIFFERENT    CONVECTION    COEFFICIENTS    WHEN   AVERAGE   ALUMINUM 
PLATE     TEMPERATURES     ARE    EQUAL 
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fact  that  the  high  conductivity  allowed  this  heat  to  be  rapidly  transferred 
to  the  surface.   It  was  also  found  that  the  temperature  of  both  a  one-half 
and  one  inch  aluminum  plate  were  nearly  identical.   Consequently,  it  was 
expected  that  the  division  of  the  plates  would  also  be  nearly  identical  since 
the  resistance  to  the  heat  flow  does  not  vary  be  a  great  amount. 

When  convection  heat  transfer  is  made  available,  the  amount  of  heat 
flow  through  the  inner  surface  of  the  outside  plate  apparently  exceeds  that 
through  the  radiating  surface  of  the  plate  by  some  considerable  amount,  ex- 
cept for  low  convection  coefficients,  as  can  easily  be  seen  in  Fig.  11  . 
Thus  it  is  entirely  expected  that  the  percentage  division  of  the  aluminum 
plate  will  remain  relatively  constant  until  the  convection  coefficient  at- 
tains a  low  value.  When  the  heat  has  been  rapidly  transferred  to  the  sur- 
face, the  presence  of  an  expeditious  heat  transfer  medium,  through  the  avail- 
ability of  convection,  allows  this  heat  to  be  removed  from  the  system. 

With  regard  to  the  maximum  thermal  stress  of  a  plate  which  is  held  in 
perfect  restraint,  Fig.  5  may  again  be  utilized.   The  temperature  at  which 
the  maximum  thermal  stress  first  exceeded  the  yield  stress  for  7075-T6  alu- 
minum was  about  790°R.   Examination  of  Fig.  11  shows  that  when  convection 
is  made  available  to  the  aluminum,  the  steady  state  temperature  does  not 
exceed  790°R  until  the  convection  coefficient  has  reached  approximately 
5  Btu/hr,ft  ,°F,  with  the  exact  value  varying  slightly  with  the  plate  thick- 
ness.  Since  there  is  a  good  probability  that  the  temperature  of  a  gas  flow- 
ing in  a  channel  such  as  has  been  postulated  here,  will  be  lower  than  the 
assumed  value  of  500°R»  and  since  the  convection  coefficient  will  almost 
surely  be  greater  than  5  Btu/hr, ft2,°F,  the  aluminum  can  be  subjected  to  a 
much  higher  flux  level  than  was  found  for  a  plate  with  radiation  heat  loss 
only,  without  endangering  its  structural  stability. 


*>2 


10  20 

CONVECTION     COEFFICIENT,  h,    (BTU/hr,  Ft.2  °F) 

FIG.  14.      RELATION    OF    MAXIMUM  THERMAL  STRESS  FOR 
STEEL    PLATE    RESTRAINED    IN    BENDING  TO  CONVECTION 

COEFFICIENT 
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Figure  1*4-  shows  the  relation  of  the  maximum  thermal  stress  in  a  steel 
plate  restrained  in  bending  to  the  convection  coefficient.   This  curve  is 
seen  to  exhibit  a  distinct  minimum  as  the  thickness  of  the  plate  increases. 
In  the  present  instance,  the  stresses  for  a  two  inch  plate  were  not  plotted. 
For  this  case  the  temperatures  for  the  lower  convection  coefficients  reached 
values  for  which  the  elastic  properties  of  steel  are  no  longer  considered  to 
be  completely  valid  (3).   The  presence  of  this  minimum  is  the  result  of  the 
availability  of  a  mode  of  heat  transfer  from  both  surfaces  of  the  plate  and 
also  the  relative  importance  of  each  mode. 

The  maximum  thermal  stress  for  a  plate  restrained  in  bending  is  found 
at  the  edge  of  the  plate  for  each  of  the  cases  considered.   Furthermore, 
this  stress  is  always  found  to  be  one  of  tension.   For  a  steel  plate  at  a 
relatively  high  convection  coefficient,  the  greater  amount  of  heat  is  lost 
from  the  inner  face  of  the  second  plate.   Consequently,  the  temperature  at 
this  face  is  lower  than  the  temperature  at  the  face  losing  heat  by  radiation, 
and  the  quantity  (T_._  -  T)  is  greater  at  the  inner  face.   Thus  there  is  also 
a  higher  thermal  stress  at  this  face.   However,  the  heat  lost  from  the  outer 
face  by  radiation  soon  begins  to  become  more  significant  as  the  convective 
heat  transfer  coefficient  is  lowered.   Thus  the  difference  between  the  temp- 
erature at  the  inside  face  and  the  average  temperature  of  the  plate  becomes 
less,  while  the  difference  at  the  outside  face  becomes  greater.   This  also 
causes  a  lowering  of  the  thermal  stress  at  the  inside  face  of  the  plate  while 
the  stress  at  the  outside  face  is  rising. 

When  the  equivalent  gray  body  convective  coefficient  is  equal  to  the 
convection  coefficient  at  the  inside  face,  the  temperatures  at  both  the  in- 
side and  the  outside  face  are  equal.   At  this  point  the  maximum  thermal 
stress  for  a  plate  restrained  in  bending  reaches  a  minimum.   For  any  lower 
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convection  coefficient,  the  temperature  difference  (Tav  -  T)  becomes  greater 
at  the  outside  face  and  continues  to  rise  at  a  rapid  pace  as  the  radiation 
heat  loss  takes  on  continuously  greater  significance.   Consequently,  the 
maximum  thermal  stress  also  rises  rapidly  as  is  shown  in  Fig.  1^. 

For  an  aluminum  plate,  Fig.  15  does  not  seem  to  exhibit  the  same  prop- 
erty just  discussed  for  steel.  However,  further  examination  of  the  equiva- 
lent gray  body  convective  coefficient  shows  it  to  be  below  a  value  of  one 
Btu/hr,ft  ,°F  for  all  the  steady  state  cases  considered.   Thus  the  minimum 
point  has  not  been  reached  within  the  range  of  convection  coefficients  in- 
vestigated. 

As  was  found  to  be  the  case  for  the  single  plate,  the  maximum  thermal 
stress  for  a  plate  restrained  in  bending  is  not  of  such  a  magnitude  that  it 
becomes  significant  in  the  structural  considerations  for  any  of  the  cases 
which  have  been  presently  examined.  Although  this  is  true,  the  presence  of 
the  minimum  which  has  been  discussed  above  should  be  kept  in  mind  in  consid- 
ering the  presence  of  thermal  stresses  in  other  materials  not  examined  in 
this  study. 

Conclusions 

When  radiation  heat  transfer  is  available  to  a  flat  plate  subjected  to 
a  high  intensity  gamma  photon  flux  with  a  magnitude  of  9^500  Btu/hr,ft  , 
the  temperatures  remain  in  a  range  where  its  structural  utility  is  not  seri- 
ously impaired.   In  the  case  of  steel,  a  temperature  of  about  2000°R  is 
reached  in  approximately  1  5  minutes  for  each  of  the  thicknesses  considered 
except  for  a  1  jk   inch  plate .   In  that  instance  the  steady  state  temperature 
does  not  reach  2000°R.   In  contrast  to  the  nature  of  this  temperature  rise, 
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CONVECTION     COEFFICIENT,    h,    (BTU/hr,  Ft.toF) 
FIG.  15.       RELATION    OF    MAXIMUM   THERMAL  STRESS    FOR 
ALUMINUM    PLATE     RESTRAINED    IN  BENDING  TO  CONVECTION 

COEFFICIENT 
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aluminum  is  found  to  approach  the  steady  state  much  more  slowly.  Hence, 
even  if  aluminum  were  exposed  to  the  same  intensity  flux  as  steel,  the  dan- 
ger of  nearing  the  melting  point  would  be  negligible  for  a  period  of  reactor 
operation  in  the  range  of  1 5  to  30  minutes. 

After  shutdown,  the  temperature  of  the  plate  cools  quite  rapidly  despite 
the  presence  of  a  time  dependent  gamma  photon  source  from  the  fission  pro- 
ducts.  The  difference  between  this  cooling  and  that  for  no  heat  generation 
after  shutdown  remains  nearly  the  same  for  a  period  of  several  minutes.   The 
continuous  nature  of  the  temperature  as  a  function  of  time  immediately  after 
shutdown  is  difficult  to  exhibit  and  is  probably  more  evident  in  the  consid- 
eration of  extremely  small  time  increments.   However,  this  is  of  little  or 
no  consequence  to  a  macroscopic  investigation  of  the  present  nature. 

When  additional  modes  of  heat  transfer  are  made  available  to  a  flat 
plate,  the  steady  state  temperature  can  be  significantly  lowered.  For  steel, 
the  additional  convection  heat  transfer  becomes  significant  at  about  1 0 
Btu/hr,ft  , °F.   At  this  point,  approximately  the  same  amount  of  heat  is  lost 
by  both  convection  and  radiation.   For  convection  coefficients  below  this 
value,  the  steady  state  temperature  rises  rapidly.   In  aluminum,  the  critical 
value  for  the  convection  coefficient  is  below  one  Btu/hr.ft  , °F.   The  rise 
of  the  steady  state  temperature  for  convection  coefficients  below  this  criti- 
cal value  is  also  considerably  more  severe  than  for  steel. 

The  only  apparent  thermal  stress  damage  which  might  occur  for  any  of 
the  cases  considered,  is  for  a  steel  plate  which  is  perfectly  restrained. 
This  physical  situation  is  somewhat  unlikely,  however,  since  almost  all  the 
structural  members  in  a  nuclear  rocket  would  be  subjected  to  a  high  inten- 
sity gamma  photon  flux  and  would  have  a  tendency  to  expand  simultaneously. 
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Although  such  a  system  would  still  remain  restrained  in  bending,  the  maximum 
thermal  stresses  of  this  type  are  small  compared  to  the  elastic  limits  of 
the  systems  examined.   It  is  significant  that  subjection  of  an  aluminum  plate 
to  a  maximum  thermal  stress  of  the  type  encountered  in  perfect  restraint, 
does  not  present  difficulties  until  the  plate  has  been  irradiated  by  the 
gamma  flux  for  a  period  in  excess  of  the  15  to  30  minutes  which  might  be 
considered  normal  for  the  operating  time  of  a  nuclear  rocket.   Furthermore, 
by  providing  even  minimal  additional  modes  of  heat  transfer,  the  steady  state 
maximum  thermal  stress  of  an  aluminum  plate  can  be  kept  well  below  the  yield 
stress. 

Finally  the  case  of  a  plate  with  heat  transfer  available  at  both  faces 
exhibits  a  distinct  minimum  for  the  maximum  thermal  stress  of  a  plate  re- 
strained in  bending.   For  all  the  metals  and  conditions  which  have  been  ex- 
amined in  this  study,  the  magnitude  of  this  type  of  thermal  stress  is  small 
enough  not  to  be  of  danger  to  structural  members.   However,  there  are  un- 
doubtedly some  metals  for  x^hich  the  minimum  exhibited  by  this  curve  could 
be  of  primary  importance  in  the  design  of  a  nuclear  rocket. 

Suggestions  for  Further  Work 

The  availability  of  further  data  in  the  areas  discussed  in  this  study 
lies  principally  in  refinements  of  calculational  techniques  and  in  the  vari- 
ation of  certain  parameters  to  obtain  a  more  accurate  description  of  the 
systems  under  consideration.   All  the  calculations  which  have  been  present- 
ed, used  the  basic  assumption  that  there  was  no  Compton  scattering  in  the 
plates  examined.   Of  course,  even  for  the  thinnest  plate,  this  assumption 
is  not  precisely  accurate.  A  more  exact  evaluation  of  the  mode  of  gamma 
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photon  absorption  could  be  obtained  by  use  of  the  moments  method  (8).   This 
could  serve  to  set  limits  on  the  thickness  of  a  plate  which  is  reasonably 
approximated  by  the  method  employed.   Also,  it  would  give  a  more  accurate 
description  of  the  temperature  response  of  materials  of  greater  thickness 
than  have  been  considered. 

Another  interesting  study  would  be  a  closer  investigation  of  some  of 
the  parameters  involved  in  the  calculations.   Rather  than  consider  some  spe- 
cific metal,  it  would  be  possible  to  choose  some  reasonable  set  of  physical 
constants  and  then  vary  each  one  individually  while  holding  the  others  steady, 
In  this  respect  it  would  also  be  possible  to  vary  the  temperature  dependence 
of  each  of  these  physical  constants  after  assuming  several  initial  values. 
Those  constants  which  would  be  of  special  interest  in  both  the  temperature 
and  thermal  stress  equations  are  the  density,  the  heat  capacity,  the  emis- 
sivity,  the  coefficient  of  thermal  expansion  and  the  thermal  conductivity. 

The  thermal  conductivity  would  be  of  special  interest  because  of  its 
effect  on  the  use  of  the  "thin  body"  approximation.  A  more  thorough  investi- 
gation of  the  magnitude  of  the  conductivity,  coupled  with  the  investigation 
of  the  thickness  of  the  plate  which  would  allow  the  use  of  the  "thin  body" 
approximation  would  prove  very  helpful.   Since  the  "thin  body"  equation  is 
an  ordinary  differential  equation,  the  resultant  temperature  expressions  can 
be  evaluated  much  more  rapidly. 

It  would  also  prove  interesting  to  define  the  physical  system  more  ac- 
curately to  determine  the  effect  of  several  other  parameters.  A  system  could 
probably  be  constructed  which  would  present  a  greater  surface  area  to  the 
environment  to  make  the  radiative  heat  loss  more  efficient.   In  the  present 
considerations,  the  view  factor  has  been  assumed  to  be  unity  in  all  cases. 
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In  varying  the  spatial  characteristics  of  the  system,  the  view  factor  would 
also  require  further  investigation.  These  alterations  should  also  prove 
helpful  in  attempting  to  remove  excess  heat  by  the  use  of  convection,  such 
as  was  attempted  through  the  use  of  two  parallel  plates  with  a  flow  of  gas 
between  them. 

Similarly,  it  would  also  be  possible  to  undertake  a  more  thorough  in- 
vestigation of  the  double  plate  case.   The  temperature  of  the  gas  flowing 
through  the  channel  between  the  plates  should  be  followed,  as  well  as  the 
variation  of  the  heat  transfer  coefficient.   In  the  same  vein,  the  variation 
of  the  sink  temperature  on  the  amount  of  heat  removed  as  well  as  the  temp- 
erature distribution  should  prove  of  interest.  This  would  be  true  for  the 
double  plate  case  as  well  as  each  of  the  other  physical  situations  which 
have  been  examined. 

It  should  also  prove  of  considerable  interest  to  more  closely  examine 
the  flux  levels  on  the  temperature  and  thermal  stresses  after  certain  pre- 
determined times  of  reactor  operation.  Both  the  magnitude  of  the  energy 
flux  and  the  distribution  of  this  energy  would  have  a  considerable  effect 
on  the  temperature  response  of  any  of  the  systems  which  have  been  considered 
in  this  study  or  those  which  have  been  discussed  above. 

Another  possibility  which  should  prove  useful  in  the  structural  consid- 
erations of  nuclear  rockets  is  the  incidence  of  a  plane  flux  beam  on  the 
side  of  a  metal  rod  or  other  cylindrical  supporting  member.  The  solution 
of  the  resultant  partial  differential  equation  in  cylindrical  geometry  and 
the  use  of  any  of  several  ordinary  boundary  conditions,  presents  a  two  di- 
mensional system  which  is  not  easily  solved  by  conventional  methods.  The 
use  of  Laplace  transforms  may  prove  useful  in  attempting  to  express  the 
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discontinuous  boundary  conditions  which  might  commonly  be  used.  However, 
the  use  of  numerical  techniques,  again  coupled  with  the  aid  of  a  high  speed 
digital  computer  would  probably  be  found  most  helpful. 

The  magnitude  of  the  maximum  thermal  stresses  found  for  the  metals  in- 
vestigated are  well  within  the  elastic  limits.   The  metals  examined  are  quite 
common,  however,  and  their  properties  should  be  subjected  to  more  thorough 
scrutiny.   An  examination  of  the  inelastic  qualities  of  these  and  other  com- 
mon metals  would  seem  in  order.   Of  special  interest  would  be  the  creep  of 
structural  members  caused  by  a  cyclic  temperature  change  which  would  occur 
when  a  reactor  is  operated  for  a  short  time  and  then  shut  down,  periodically. 
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Derivation  of  the  Temperature  Distribution  and  Thermal 
Stresses  in  a  Flat  Plate  with  Insulation  at  One  Face 
and  Radiation  Heat  Loss  at  the  Other  Face 


The  partial  differential 
equation  for  the  temperature  dis- 
tribution in  the  flat  plate  shown 
in  Fig.  16  is 


£&    .H=j_<2e 


a**  ^*  "t>f at 


(A-1) 


*-x 


FIG.  16. 
CO-ORDINATE    SYSTEM 


The  heat  generation  term,   H,    is 
given  by  Bonilla   (3)   as  yule->ttx, 
where   the  absorption  coefficient, 
JU,    is   such   that   it  considers   +he  effect   of  buildup   through   the  plate.      Al- 
though  this  is  not  truly  a  constant,    it   can  be  represented   in   that  manner 
for   thin  systems   such   as  are  being  considered  here.      The  boundary  conditions 
for  this  system  are 


-*«L-» 


--»ISL^»L 


and  the  initial  condition  is 


(A-2) 


(A-3) 


ei? 


(A-*0 


where  hr  is  the  equivalent  gray  body  convective  radiation  coefficient, 
defined  by 
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-£=k(T-T.)  =«r(T'-T<)  (a-5) 

The  assumption  is  first  made  that  6  may  be  written  as 
Then  operating  on  Eq.(A-6)  with  (A_1 )  yields 


?*«   +  a**         ^           of  3t  lA~7; 

Epel  (5)  shows  that  such  an  equation  may  be  solved  by  separating  this  equa- 
tion into  two  equations  in  X  and  X"  alone  and  setting  them  both  equal  to 
zero.   This  leaves 

JSLJL   +£-l£*l     =0  (A.8) 

d  %x         Jt 

Integrating  Eq.(A-8)   twice,    the  solution  becomes 

x  =  -5jr-  ■|"c'%  +  c*  (a-1o) 

In  order  to  solve  for  the  constants  C-|  and  C?  it  is  necessary  to  transform 

the  boundary  conditions  so  they  are  compatible  with  the  system  of  equations 
being  solved.   For  Bq.(A_8)  these  are 

-A^r)    =0  CA.11) 

->#)  =kX)  GUI*) 
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In  the  case  of  Bq.(A-9)  they  become 

-A$EL-°  (A-13) 

The  initial  condition  now  becomes 

T  =  -X  (*-i3) 

Applying  Bq.(A_1l)  and  (A-1 2)  to  (A-10),  there  results 

C,  ="^-  (A-16) 

C*  -  JL(\-e^L)  +i.(^L+e^L)  (1.17) 


Thus 


(A-18) 


which  is  the  steady  state  temperature  distribution.   Bq.(A-9)  may  be  solved 
by  again  separating  variables  in  such  a  manner  that 

TaU(t)Sw  (A-19) 

Operating  on  Bq.(A-19)  with  (A-9)  yields 

13^.  =   _!_^LL  (A-20) 

■     5  3p?        3iU3t 

This  is  a  familiar  eigenvalue  problem  and  has  solutions  only  when  both  sides 
of  the  equation  are  equal  to  the  same  constant  —  X   .   These  solutions  are 
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U=e"  (A-21) 


5  =  C3s^Xx    +C,Co-,A*  (A-22) 

From  which 

/  i  \        *Xt 

'X  -  (C3^*a-*    v  C4  c^-,A*)e."  (A-23) 

Applying  Bq.(A-13)  to  (A-23),  Cj   is  found  to  be  zero.   Then  applying  Bq. 
(A-14)  to  (A-23),  the  eigenvalue    is  given  by 

There  are  a  countable  infinity  of  roots  (  KJ   to  this  transcendental  equation. 
Using  this  information  and  applying  Bq.(A-15)  to  (A-23)  gives 

-£  C«(~*  X.*X,)  =  -^(e'1-  e-**)  +  £(l_-  *)+  i-( \  -  e*)      (a-25) 

Multiplying  both  sides  of  Eq.(A-25)  by  the  orthogonal  function  c-o-s  X„* 
and  integrating  over  the  period  of  orthogonality  yields,  after  some  trigo- 
nometric manipulation 


From  which 


(A-26) 


7^__  ZI/cj.6^  [VU*A»  4-Q   +  A*(H.*  -^eT^gcA-L 

L      ^Ai(Alnv)  -*zA!l  +w;l  -*-m< 


and   the   temperature  is   found   to  be 


e 


(A- 2?) 
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1   *At 


(A- 28) 


which  satisfies  the  assumption  concerning  the  separability  of  Q. 

It  is  possible  to  arrive  at  this  same  equation  by  transforming  the  par- 
tial  differential  equation  and  the  boundary  conditions  from  the  time  domain 
to  the  Laplace  domain.   Bq.(A-l)  then  becomes 


ft  -^^  =^*s  -M 


Q**      'Its  7{ 

where  8  is  the  transform  of  S  and  also  where  0(0)  is  zero. 
The  transformed  boundary  conditions  are 

Applying  these  conditions  results  in  the  equation 


(A- 29) 


(A- 30) 


(A-3D 


h„t*»V(Ll55t)  -^(Lj^)^i^k(LV^y 


(A-32) 


which,  when  inverted  back  into  the  time  domain,  yields  Bq.(A_28)  above. 

Timoshenko  and  Goodier  (21)  show  that  for  a  temperature  distribution 
in  the  x  direction,  the  thermal  stresses  in  the  y_  and  z   directions  are  given 
by 
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?  =  <*T  -  fl-f"  9  +  L^  +  'z  faJ  ^(y  -  Vfe)  B<A*}        (A-33) 

The  first  term  within  the  brackets,  <j-    ,    is   simply  the  negative  of  Eq.(A-28) 


The  second   term  becomes 


_  r 


^"E 


*/ 


J  *  Wrl  " 


where  Gn  is  defined  by 


>  (A  A\  +  C)  +  XI  ( M  ^le^ec  U 


y?)  JKX\\L  +k:l  +-M* 


] 


(A-35) 


Similarly,  the  third  term  is 


jr*^  fee  -")*>(="-  o 


rts»     i- 


-XL* 


(A-36) 


^>v<w 


AnL 


Combining  these  three  terms  yields  the  final  equation 


^  .^ob**#*^5Cr  ♦  *^-  ^ 


(A- 37) 
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APPENDIX  B 


Derivation  of  the  Temperature  Distribution  and  Thermal 
Stresses  in  a  Flat  Plate  with  Decaying  Internal  Heat 
Generation 


The  partial  differential  equation  describing  the  conditions  for  this 
physical  case  is 

a**         jt        ~o{  at  (B-1) 

where   the   time  dependent    flux,    l(t),    is  given  by  Lottes    (12)   as 

I(t)  =  I(o)  jo.i(t  +  »o)  °*oo8T(t  ♦  *-  lo'f  -  o.oots  e" (     -  o.oo  i^e  C*==J 

,    *    v  (B-2) 

ol(t^t%+-)o/>  -o.o5T(t^t4+2Mo7  -ooozse        -  0.0013 <s  v?9<»~oy     | 

where  t  is  the  time  after  shutdown  and  t_  is  the  time  of  reactor  operation, 
both  times  being  in  seconds.  In  such  a  form  Eq.(B-l)  is  not  readily  solved 
by  the  separation  of  variables  technique.  If  the  Laplace  transformation  is 
utilized,  the  resulting  equation  in  the  Laplace  domain  is 


.*•=• 


§£  .  £(»6  -  b<«>)  -£  e-*X  { ita]  (b-3) 

where  0  is  the  transform  of  S.   The  initial  condition  B(0)  is  given  by 

and  *T{l(t)j   is  the  Laplace  transform  of  Eq.(B-2).   The  transformed  bound- 
ary conditions  are 


-**  5xL0= 


(B-5) 
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-j$L  =K$~  (B-6) 

where  hr  is  again  the  equivalent  convective  radiation  coefficient,  given  by 


l=k(T-T.)=^r(T«-T„<) 

The  solution  to  this  system  of  equations  is 


(B-7) 


•  -  coKClJ^)  -  *l£5£  ^(0*^)  1 


(*-V 


JL 


jtjv*  •sao^u***;  -^  ^  cosv  u.\Va) 


h,  Cos  a-*l.  w",l 


C"7    Ci,0?.  W„  cosA«L     ~AX»SlA<.X>L 

■UPS* 


to 


-sAx^f^l) 


where  Gn  is  given  by  Bq.(A_35)  and  also  where 

~f  $Tl+\l   -Ti  cf*°*i*'c:  l*l\       QQgTfe^o1)'"  Q.QOZ5      __  QQO)3 


4- 


+ 


(B-8) 


ooo 

(B-9) 


-OX 


^z.o4o 


S  -»-    |/z.o0# 


ooo 
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in  which  s1    is  the  Laplacian  variable  corresponding  to  the  time  t1  =  t/1 0  and 
sw  is   that  corresponding  to  the  time  t*  =  t/(tg  +  10).     In  addition 

E*(s)=ye"ff"*^f  (B.10) 

• 
It  is  now  necessary  to  take  the  inverse  of  Eq.(B-8)  in  order  to  transfer  the 
equation  back  into  the  time  domain  and  find  the  temperature  distribution. 
Although  this  is  straight  forward  theoretically,  practically  the  fractional 
orders  in  Bq.(B-9)  present  difficulty.   Therefore,  a  satisfactory  solution 
was  attempted  by  resorting  to  the  use  of  numerical  techniques. 

Richtmyer  (18)  presents  a  discussion  of  the  methods  available  and  their 
limitations.  Written  in  finite  difference  form,  Eq.(B-l)  becomes 

(B-11) 


At 


where  14J  is  a  real  constant  0  -  V  -  1 .   For  a  known  temperature  distribution 
at  time  nAt,  this  system  may  be  solved  for  the  (n  +  1  )st  time  increment  at 
any  given  x  =  jax  for  each  of  J  increments.  The  constant  ¥   is  chosen  as 
the  value  which  eliminates  first  order  errors  and  still  retains  absolute  sta- 
bility. This  value  is 

The  boundary  condition  at  the  inside  face  is 

-*$)     -o  CB-13) 


7^ 


which  would  be  written  in  finite  difference  form  as 


-X 


8.  ~  Bo 


£% 


=  o 


(B-1i*) 


However,  Eq.(B-13)  is  really  a  symmetry  boundary  condition  and  it  is  more 


correct  to  write  this  equation  in  finite  difference  form  as 
r  -A 

-A 


s,  -  a, 


Ay. 


The  other  spatial  boundary  condition  is 


(B-15) 


-k 


ft  -  ar.. 


=  K% 


(B.16) 


The  initial  condition  is 


b;  =  8(^ 


(B-17) 


at  the  time  of  shutdown.  Since  the  distribution  at  time  nAt  is  known,  it 
is  possible  to  use  Bq.(B-H)  to  set  up  a  series  of  equations  which  may  be 
solved  for  the  temperature  distribution  at  time  (n  +  l)At.  Using  ten  in- 
crements, this  may  be  put  in  matrix  form  as  written  symbolically  below 


~2{aF 


+  2.4J 


+  ?_*+» 


b3    -*v&~:   =17     i»i 


^#  ^H^(r™)l  -rap  -13  *«» 


(B-18) 
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It  is.  also  necessary  to  express  the  thermal  stress  distribution  in  a 
free  plate  ><y  means  of  some  numerical   technique,    since  an  analytical  expres- 
sion for   'ne  temperature  distribution  is  not  available.     This  may  be  accom- 
plished by  expanding  Bq.(9)   through  the  use  of  Simpson's  rule.      Thus,    it  is 
found   Jiat  when  the  plate  is  divided  into  ten  equal  increments,    there  results 

^U  =(a+  4Bt+*.8JH&,+  •••  +  *6,  +  48L+8.)(fe)    (B.19) 

o 
/nd 

(B-20) 

whence,    the  final  expression  for  the  thermal   stress  in  a  free  plate  at  each 
of  the  j  mesh  points  becomes 

<^  =  <S  =-8-"+55-(&,  +  48,+ze,+484  +2es+4e1+za,-M&8+2V4Qi(*e.,) 


+° 


(B-21 ) 
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APPENDIX  C 


Derivation  of  the  Temperature  Distribution  in  Parallel 
Plates  with  a  Constant  Gamma  Photon  Flux,  Insulation  at 
the  Inner  Face,  Gas  Flow  Between  the  Plates  and  Radiation 
Heat  Loss  from  the  Outside  Face 


The  inside  plate  presents  essentially  the  same  problem  as  was  solved 
in  Appendix  A.  However,  in  this  case  the  heat  loss  from  the  outer  face  of 
the  plate  is  by  convection  to  a  gas  flowing  between  the  two  plates.   This 
boundary  condition  is  then 

-X^r)     =K6')  L  (c-1) 

where  the  subscript  one  is  used  to  denote  those  quantities  particular  to 
the  first  plate.  Similarly,  the  subscript  two  will  be  used  for  the  second 
plate. 

The  solution  of  Bq.(A-l)  using  the  boundary  conditions  given  by  Eq.(A-6) 
and  (C-1 ),  and  the  initial  condition  given  by  Eq.(A-10)  is 

(C-2) 


IP 


n«i 


A  U,  X^,  ±   k)  *  \n.  (  ^  A,  ~A]m)^  4ec\W,L. 


In  the  case  of  the  second  plate,  the  flux  incident  on  the  inside  face 
is  given  by 

I,  =  I,  e-^  (c-3) 

The  temperature  distribution  is  again  expressed  by  the  partial  differential 
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equation 

The  assumption  is  also  again  made  that  this  equation  can  be  separated  into 
the  form 

6,=  XW  +l>.,t)  (c-5) 

Operating  on  Bq.(C-5)  with  (C-*0  and  then  separating  the  resultant  into 
two  equations  in  Xp  and  tp  an<^  setting  each  of  these  equal  to  zero  again 
leads  to 

These  equations  are  solved  in  the  same  manner  as  was  shown  in  Appendix  A, 
and  yield  the  equations 

Tx  «  (c-'Vuvl^*-*  +Cmt**Atx,Je  (c.8) 

X^  ~  X\e"^^    +  £3**    +  ^4  (c-9) 

The  boundary  conditions  for  the  second  plate  are 

-*MA      =  kB)  (do) 

where  h  is  again  given  by  fii.(A-13). 
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'0 

which  now  becomes 


For  Eq.(C-8)  the  transformed  boundary  conditions  are 

-^^ri^^U  (c"13) 

For  Bq.(C_9)  they  are- 
The  initial  condition  is 

't«0 


(C-16) 


Applying  Eq.  (0-1*0  to  (C-9)  gives 

Mi+C>)-K^+C<)    ■  (c-18) 

And  applying  Bq.(C-15)  to  (C-9)  gives 

^(^^^)=K(-^-hC,u  +  c4)      (,19) 

These  last  two  equations  may  be  solved  simultaneously  to  yield 

5  JkA    +>,h,    *KHTU 
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^-^Hfc.+*fc) 


MvJsJiZf!!^  ■H.Oi.eT'^  +  h.V 


This  reveals  the  steady  state  solution  to  be 


+  V-* 


(C-21 ) 


*w^0-e-^)  -^(ke^+yj 


Then  substituting  Eq.(C-8)  into  (C-12)  gives 

1   X,  ^ 


U-*,+^l 


(C-22) 


(C-23) 


And  operating  on  Eq.(C_8)  with  (C-13)  and  making  the  substitution  establish, 
ed  in  Eq.(C-23)  there  results  the  transcendental  equation 


(C-2*0 


where  Xjjo  represents  the  countable  infinity  of  roots  to  this  equation. 
Carslaw  and  Jaeger  (*0  show  that  if  the  steady  state  solution,  X  =  fCxg), 
can  be  developed  in  an  infinite  series  such  that 


too  =£c»JL 


(C-25) 


where  X.  i*  °*"  th®  same  form  as  the  spatially  dependent  part  of  t*2»    then 
the  solution  to  t*  2  is  given  by 

_*,X.*t 


(C-26) 
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where  C^  is  found  from  the  orthogonal  properties,    and  is  given  by 


(C-27) 


0  -    -HX 

and  the  denominator  remains  to  be  evaluated. 


From  the  solution  for  the  spatially  dependent  part  of  "C^ 


■»*. 


(C-28) 


This  is  first  put  into  the  form 


x. fx>.  -  -f X..  ^ *. 


(C-29) 


from  which  this  is  evaluated  as 


x.fx>.«- 


X. 


«>.to'- 


(C-30) 


Substituting  Bq.(C_23)  into  the  spatially  dependent  part  of  (C-8)  and  mul- 
tiplying both  sides  of  the  equation  by  A^*  results  in 


Taking  the  derivative  of  X^  yields 

Squaring  the  above  two  equations,  there  results  immediately 

Operating  on  both  sides  of  this  equation  by   J^  dLy,  gives 


(C-31  ) 


(C-32) 


(C-33) 


MV^fGHM^feJ] 


L. 


(0-3*0 
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Adding  B*.(CL30)  to  (C-3*0  produces 

Bq.(C-14)  may  be  rearranged  to 
dX,v  _  Ju  V 


*»«*•» 


n,«o 


d^t      *T 


and  then 


V*   _* 


(ft?  -ftjx: 

which  shows,    from  E*l.(C-33)»    that 


x:=i 


and  also  by  substitution  from  (C-36),  at  Xp  =  0 


v    dX»» 


-* 


Similarly,    Bq.(C_15)  becomes 


Ju 


»M 


from  which 


&■)'  =(t)'x: 


X. 


(C-35) 


(C-36) 


(C-37) 


(C-38) 


(C-39) 


(CJ*0) 


(C-41) 


Substituting  this  into  Eq.(C_33)  shows  that 

Then  substituting  first  from  Bcl.(CJK))  and  then  from  (C-^-2),  there  results 


(CJ*2) 


at  xp  =  Lp 
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MH =-&)(****-) 


(CJO) 


Consequently, 

AX. 


X. 


eLjtj 


'.^.   r 


t 


(C-^4) 


Then  from  Eq.(C-35),  there  results 


X>.  -4j-it(W  ♦"*  gfetflft1*  t 

(CL*5) 

For  the  evaluation  of  the  numerator  of  Eq.(C-27),  Bq.(C-22)  is  rearranged 
to  the  form 


too  =  xz = n****  ->-"Rm^  -m. 


(CJ*6) 


where 


*-& 


K-— 


i. 


(c-47) 


Performing  the  indicated  integration  on  Bq.(C-27)  yields 
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W-^  cos  X^.  +  *i~ A„^)  (*R.  e*"*  +-  1^-**  -v-  "R J  oUz 


=  ^-  j(co*X^  +  X^*  si^X^J  -  j^-  6^**^0*  A^%*-  s</nX^ 


*i*L» 


»»•• 


+T4       jx~^) 


(C-48) 


+ 


[_eZ^l 


(•^g^sU^X^t"^t  -  Anx  cos  X »»•*•) 


(^  +/) 


■]} 


*.^» 


*te« 


*1^  (~*X  ^^-^  -f  sUn.X.,,0 

Which  may  be  separated  and  written  as 

a 


%i=L, 


ht.«o 


J 


"«*. 


-BS^fc) 


+  SUM. 


un.  \„%LX 


[x.(^^)-5fe.-fe$,(x..-ffc)]<«. 


JttX\.       XV       (£~v?J       >.  U*.  +/.) 


Then  substituting  (0-27)  into  (0-26)  to  solve  for  T2,  there  results 
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6t= -R.e--"""  +-R.*.  +-R3 


+SF^[-tt.+%.(,-iti-ajS^*» 


+-  ^[t^(-t-&)+s^^-^a] 


(C-50) 


♦^♦Ht-^Hie*-^ 


+  Oos  X„/¥ 


} 


\ 


*x. 


.}'] 


Considering  the  case  where  h£  =  0,  this  solution  should  reduce  to  the 
solution  obtained  in  Appendix  A.  This  is  found  to  be  the  case  and  is,  in 
part,  a  verification  of  the  present  solution. 
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APPENDIX  D 


Description  and  Explanation  of  the  IBM  1620  FORTRAN 
Program  to  Find  the  Temperature  Distribution  and  the 
Thermal  Stresses  in  a  Plate  with  Insulation  at  One  Face 
and  Radiative  Heat  Loss  at  the  Other  Face  when  Subjected 
to  a  Constant  Gamma  Photon  Flux  at  the  First  Face 


The  equation  for  the  temperature  distribution  which  must  be  solved  is 

(D-1) 

L  *xux>y)L rx\L^iL^jkK F 

n*i 

In  order  to  find  a  numerical  solution  to  this  equation,   it  is  first  necessary 

to  complete  two  iterations.     The  first  is  for    An  which  is  found  from  the 
transcendental  equation 

AAa/nAj-  ^jj*-  (D-2) 

It  is  also  necessary  to  find  the  correct  value  for  hr,  which  is  a  function 
of  the  surface  temperature  and  is  given  by 


W_fe<r(T4-T.4) 


as  well  as  by  Eq.(D-l)  above. 

It  is  possible  to  carry  out  these  two  iterations  simultaneously.  A 
value  for  hr  is  first  assumed,  from  which  each  of  the  X*s  may  be  found. 
For  the  physical  cases  considered  here,  it  was  found  that  the  first  two  solu- 
tions to  Eq.(D-2)  were  sufficient  to  give  accurate  answers.  Using  these  val- 
ues, 9  may  then  be  calculated  from  Eq.  (D-1  ) .  The  temperature  obtained  is 
used  in  Eq.(D-3)»  and  this  calculated  value  for  hr  is  compared  with  the 
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assumed  value.  Before  a  reasonable  number  can  be  chosen  for  the  next  trial, 
however,  it  is  advisable  to  further  investigate  Eq.(D-l)  and  (D-3).  Figure 
17  gives  a  qualitative  description  of  these  two  equations.  It  can  be  seen 
that  if  point  a  is  the  first  approximation,  calculation  with  Eq.(D-l)  will 
bring  the  system  to  point  b.   This  value  of  hr  can,  then  be  used  as  the  next 
estimate  from  which  a  new  temperature  is  found,  as  well  as  a  new  hr  (point 
d).   It  becomes  rather  obvious  from  Fig.  17  that  the  continuation  of  this 
process  would  take  considerable  more  time  to  find  a  reasonably  close  estimate 
for  the  true  hr  as  time  increased  and  the  temperature  rose.  However,  the 
iteration  process  is  easily  ammended  at  this  point  to  shorten  the  time  nec- 
essary to  obtain  a  satisfactory  answer.   The  graph  used  is  that  shown  in 
Fig.  18. 

The  value  for  hr  is  the  assumed  value,  and  A  hr  is  the  difference  be- 
tween this  assumed  value  and  the  value  finally  calculated  from  Bq.(D_3). 
Point  a  is  the  first  assumed' value  and  the  difference  between  a  and  b  is 
Ahr.     Point  £  is  the  second  assumed  value  and  the  second  Ahr   is  from  c  to 
d.  Point  e   denotes  the  true  value  of  hr.  A  linear  approximation  can  be  used 
to  estimate  the  point  at  which  the  true  hr  exists  ( Ahr  =  0).  This  point, 
f,  is  then  used  as  the  first  assumption  in  the  next  set  of  trials  which  pro- 
ceed in  exactly  the  same  manner.  The  convergence  to  five  significant  fig- 
ures, which  was  deemed  sufficient  for  these  calculations,  is  quite  rapid. 

When  the  best  value  of  hr  is  found,  it  is  possible  to  calculate  the  tem- 
perature distribution  and  the  thermal  stress  distribution  in  the  plate. 
Since  many  of  the  "constants"  used  in  these  equations  are  truly  temperature 
dependent,  it  is  necessary  to  adjust  these  with  each  iteration.  A  linear 
approximation  for  these  constants  was  found  (2)  and  the  average  temperature 
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of  the  plate  was  used  in  the  adjustment. 

In  calculating  the  thermal  stresses,  it  was  found  that  the  use  of  float- 
ing point  arithmetic  in  the  FORTRAN  program  caused  a  loss  of  significant  fig- 
ures when  Bq.(9)  was  used,  because  of  the  subtraction  of  almost  equal  numbers 
which  were  generated.   Consequently,  the  following  approach  was  taken  to  over- 
come  this  difficulty.  . 

The  thermal  stress  equation  for  a  free  plate  is 

t   *     \-t\Ajx  I ^e      AA       u^T'+  LguD        J 

(d-<0 

Consider  first,  the  transient  portion  of  this  expression.  The  first  two 
terms  within  the  square  brackets  can  be  expanded,  using  a  Maclaurin  series, 
to  yield 


(D-5) 


where  y  =  x/L.  Four  terms  here  were  sufficient  for  the  systems  under  consid- 
eration. The  latter  part  of  the  transient  portion  of  the  equation  is  expand- 
ed in  the  same  manner.  After  some  manipulation,  this  becomes 

Now  considering  the  steady  state  portion  of  Bq.(lX-*Ot  the  quantity  within 
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the  square  brackets  may  be  rearranged  to 
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These  last  two  terras  are  expanded  to 
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the  general  expression  for  which  is 

The  first  two  terms  within  the  square  brackets  of  Eq.(D-7)  yield 


and  the  general  expression  is 


•«»*''     V*M-», 


b(l+e^)+;l(e^,)=^)(-0   (fan-as*) 

Whence  the  final  stress  equation  is  written  as 


(D-7) 


(D-8) 


(D-9) 


(D-10) 
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(D-12) 


When  the  final  values  for  h  and  the  temperature  distribution  and  ther- 
mal stresses  have  been  determined  at  a  given  time,  it  is  possible  to  speed 
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the  iteration  for  hr  at  the  succeeding  time  by  using  the  finite  difference 
approximation  for  the  subsequent  temperature  and  then  using  the  radiation 
boundary  condition  to  calculate  hr.  The  final  equation  is 


'  u  t.-t.        +      t,-t. 


t,-t,      t,-t. 


(D-13) 
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t.*-x*    x'-t.' 


t.-t, 


— -X   .    _-_!    • 
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t,-t,   trt.J 


} 


where  the  knowledge  of  the  outside  temperature  of  the  plate  at  four  previous 
times  is  necessary  to  satisfy  the  equation.  However,  somewhat  less  accurate 
approximations  can  also  be  found  when  only  two  or  three  of  these  previous 
temperatures  are  known  by  using  either  the  first  two  or  first  three  terms 
of  this  expression. 

The  program  also  includes  several  sense  switches  which  allow  control 
of  the  output  data.  When  switch  2  is  turned  ON  during  the  time  of  the  iter- 
ation for  hr,  the  output  data  will  contain  one  line  with  the  transient  and 
steady  state  portions  of  the  free  plate  thermal  stress  respectively.  This 
is  followed  by  a  second  line  containing  the  distance  from  the  origin,  tem- 
perature, free  plate  stress,  and  steady  state  portion  of  the  tempera+ure. 
This  data  is  printed  for  each  of  the  desired,  number  of  increments  across  the 
plate.   The  last  line  contains  the  constants  X<  ,  Ao»  G-j  and  G^  respective- 
ly. If  the  switch  2  is  truned  OFF,  the  line  containing  the  transient  and 
steady  state  portions  of  the  free  plate  stress  are  omitted  and  the  free  plate 
stress  in  the  succeeding  line  is  printed  as  zero.   The  remainder  of  the 
printed  output  remains  the  same. 

During  the  running  of  the  program,  switches  1  and  3  ar©  set  in  the  OFF 
position.  When  switch  3  is  turned  ON,  the  program  will  complete  the  data 
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point  it  is  presently  in  the  process  of  calculating.  After  printing  the 
data  discussed  above,  the  constants  and  other  data  which  were  used  as  input, 
will  be  punched  in  whatever  form  they  are  presently  in  use  in  the  program. 
These  cards  are  then  available  for  use  as  the  first  several  cards  of  input 
data  when  it  is  desired  to  continue  the  program  at  a  later  time.  This  essen. 
tially  allows  continuous  running  of  the  program,  even  if  an  interruption 
should  occur.  Since  the  primary  purpose  of  these  programs  is  to  find  the 
continuous  temperature  response  of  the  plate  during  reactor  operation  and  al- 
so during  the  time  after  shutdown,  switch  1  may  be  turned  ON,  and  the  addi- 
tional' data  necessary  to  the  finite  difference  program  will  be  punched  out. 
Thus  this  program  is  both  continuous  within  itself,  and  also  completely  com- 
patible with  the  finite  difference  program  described  in  Appendix  E. 

If  either  switch  3  or  switches  1  and  3  are  ON,  the  program  will  stop 
when  these  data  cards  have  been  punched  out.  If  it  is  desired  to  obtain  the 
punched  card  output  described  above,  and  also  to  continue  the  program  after 
obtaining  these  cards,  switch  b   is  turned  ON  and  the  program  proceeds  in  the 
normal  manner.  This  switch  essentially  allows  one  to  obtain  continuous 
punched  card  output  for  any  number  of  times  at  which  the  reactor  might  be 
shut  down. 

The  input  data  cards  to  this  program  contain  the  following: 

1.  The  coefficient  of  thermal  expansion,  the  base  temperature,  an 
initial  estimate  of  the  value  0$  and  the  thermal  conductivity. 
(1  card) 

2.  The  temperature  dependent  intercept  of  the  emissivity,  the  number 
of  increments  being  considered  across  the  thickness  of  the  plate, 
an  initial  estimate  of  the  quantity  oCE/(1  -~0  )  and  the  equivalent 
gray  body  convective  radiation  coefficient.   (1  card) 
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3.  The  temperature  dependent  intercepts  of  the  modulus  of  elasticity, 

the  coefficient  of  thermal  expansion,  the  thermal  conductivity  and 

the  density  of  the  material  at  *K)°F.   (1  card) 
^.  The  temperature  dependent  slopes  of  the  modulus  of  elasticity,  the 

coefficient  of  thermal  expansion,  the  thermal  conductivity  and  the 

value  of  the  heat  capacity.   (1  card) 

5.  The  initial  values  of  the  modulus  of  elasticity,  Poisson's  ratio, 
the  time  of  operation  at  which  the  calculation  is  started  and  the 
temperature  dependent  slope  of  the  emissivity.   (1  card) 

6.  The  time  increment  between  two  successive  sets  of  calculations, 

an  estimate  of  the  value  of  h  at  steady  state  and  the  inside  tem- 
perature at  the  time  of  the  previous  calculation.   (1  card) 

7.  The  four  previous  outside  temperatures  and  the  times  at  which 
these  occurred.  The  earliest  time  considered  has  the  subscript 
one,  and  each  of  the  other  data  points  available  follows  it.   If 
less  than  four  data  points  are  available,  the  latter  data  cards  for 
which  numbers  are  not  available  can  be  set  at  any  number  and  data 
set  8  may  be  appropriately  set  to  allow  the  program  to  note  this 
fact  and  proceed  in  the  proper  manner.   (^  cards) 

8.  The  first  number  is  that  mentioned  above,  and  is  set  at  the  number 
of  previous  outside  temperatures  known.  It  is  assumed  that  one 
data  point  is  always  known  and  that  in  the  absence  of  the  correct 
temperature  this  value  is  set  at  the  base  temperature  at  zero  time. 
The  number  in  the  card  allows  the  calculation  of  the  steady  state 
temperature  distribution  before  the  time  dependent  solution  pro- 
ceeds if  it  is  set  equal  to  zero.   If  the  number  one  or  any  higher 
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number  is  inserted,  the  time  dependent  solution  proceed  immediately 
at  whatever  time  desired.   (1  card) 
9.   The  flux  and  the  absorption  coefficient  for  each  of  the  eight  ener- 
gy groups  considered.   (8  cards) 
10.   The  thickness  of  the  plate  being  examined.   (1  card) 
The  flow  chart  for  this  program  and  the  FORTRAN  source  deck  statements 
appear  on  the  following  pages. 
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1  FORMAT (El 5. 8. El 5.8 »E 15 .8.E 15.8  ) 

DIMENSION  FLUX(8) tJ(8) »VALUE(2) »ALAM(2) .OUT (4) »t(4) »GRB(2) .OLD  I  11) 
250  FORMAT!  12.12) 

read1. alpha. to .akap.cond 

readi.epsb.div.stcof.hr 

read1 .elstb.alphb.condb.densb 

read1.elstm.alphm.c0ndm.cp 

read1.elst.v.time.epsm 

read1.delt.hrstd.temp 

D0256I*1.4 
256  READl.OUTt I) .T(  I) 

READ250.M.LL 

D0255I»1»8 
255  READ1»FLUX( I ) »U( I) 

2  READ1.EL 
KK-1 
LLL»0 
|CKK«0 
EPS»EPSB 
TBSAV=TO 

SIGMA=1.7123E-09 
PI=3. 1415927 

3  TIME*TIME+DELT 
11*11 
PREV=TEMP 
MM*0 
IF(1-LL)3G0.3Q0.54 

4  N=l 
TRY1»0.1 
TRY2«1. 

5  DIFF1=HR/(C0ND*TRY1)-(SIN(TRY1*EL)/C0S(TRY1*EL) ) 
DIFF2=HR/(COND*TRY2)-(SIN( TRY2*EL) /COS ( TRY2*EL ) ) 
CURVE«(DIFF1-DIFF2)/(TRY1-TRY2) 
TRCPT=DIFF1-CURVE*TRY1 
SETM-TRCPT/CURVE) 

TEST=SET-TRY2 

IF(TEST)10.40.15 
10    TEST*(-TEST) 

15  IF( (TEST/SET)-l.E-06)40.40.20 
20  TRY1*TRY2 

TRY2«SET 

60  TO  5 
40  ALAM(N)*SET 

VALUE(N)»SET*EL 

IF(2-N)54»54.42 
42  N-N+l 

TRY1=0.1+PI/EL 

TRY2»1.+PI/EL 

GO  TO  5 
54  X*EL 

50  IF(1-LL)58»58»75 
58  DO60J=l»2 

GRB(J)*0. 
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DO601=l»8 

TRC0N=2«*U( I )*FLUX( I ) / ( COND* ( ALAM( J ) **2 )* ( ( ALAM( J )**2 ) +  ( U( I )**2) ) ) 

TRN1*U( 1 )*( (HR**2>  +  (  (ALAM( J)**2)*(COND**2  )  ) ) 

TRN2=(ALAM(J)**2)*( ( COND*HR ) - ( U ( I ) * ( COND**2 ) ) ) 

TRN3*(TRN2)*EXP( ( -U ( I ) )*EL) /COS ( VALUE ( J ) ) 

TRN«TRN1+TRN3 

TRD=HR**2*EL+HR*COND+CGND**2*EL*ALAM( J)**2 

TREXP=AKAP*( ALAM( J)**2)*TIME 

IF(180.-TREXP)61»61t60 

60  GRB(J)=GRB( J ) +TRCON*TRN*EXP I -TREXP ) / TRD 

61  SUM=0. 
D063I=1»2 

63  SUM=SUM+COS( ALAM( 1 )*X)*6RB( I ) 
70  STDY*0. 

D073I=1»8 

STD1*FLUX( I)*(EXP(-U( I )*EL )-EXP ( -U ( I ) *X ) )/(COND*U( I ) ) 

STD2=FLUX( I )*(EL-X)/COND 

STD3*FLUX( I )*(1,-EXP(-U< I ) *EL) ) /HR 

73  STDY=STDY+STD1+STD2+STD3 
IF(STDY-SUM) 110»110»77 

77  STDY=STDY+TO 

IF(  l-LL)74.74t400 

74  TEMP=STDY-SUM 
IF(EL-X)100»l00t80 

75  TPART*0.0 
D076I=1 t8 

PARTT=FLUX( I )*( 1«-EXP(-U( I  )*EL)  ) 

76  TPART=TPART+PARTT 

78  TEMP=( (TPART/(EPS*SI6MA) )  +  TO**<f ) **0.25 
GO  TO  100 

80  IFtSENSE  SWITCH  2)81t84 

81  SIGST=0. 
Y*X/EL 
BEND*Y-0.5 
D02181*l»8 
RR»1. 
SS*2. 
PP«6. 
STR2*0. 

UL«(U( I )*EL)**2 
YY=Y*Y 

R«2. 

200  S-R+l. 
P«R+2. 
RR«RR*(-R) 
SS»SS*(-S) 
pp«pp*(-P) 

STR11»6«/SS-12«/PP 

STR1«UL*(-YY/RR+BEND*STR11+1./SS> 

STR2»STR2+STR1 

SAVE-STR2 

IF(STR1)201»202»202 

201  STR1«(-STR1) 
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202  IF(SAVE)2G3.2U4#204 

203  SAVE*(-SAVE) 

204  IF(STR1/SAVE-1.E-06)215»215»213 
213  R-R+l. 

UL«UL*U( I )*EL 

YY=YY*Y 

GO  TO  200 
215  STR3»STR2*FLUX( I )/U( 1 ) 
218  SIGST»SIGST+STR3 

S1GST=SIGST*STC0F/C0ND 
220  STSUM=0. 

D0225I*1.2 

STl«(Y/2.-l./12.-Y**2/2.) 

ST2«(Y/5.-l./20.-Y**4/4.)*VALUE( 1 )**2/6. 

ST3=(Y/10.-l./36.-7.*Y**6/45. )* VALUE  I I )**4/112« 

ST4«(90.*Y**8-8.-4«*Y)*VALUE( I ) **6/3628800. 
225  STSUM=STSUM+GRB( 1 )*VALUE( I ) **2* ( ST 1-ST2+ST3+ST4 ) 

STSUM*STCOF*STSUM 

STRSS=SIGST+STSUM 

PRINT1»STSUM»SIGST 

GO  TO  86 

84  STRSS*0. 

86  PRINT1»X.TEMP»STRSS»STDY 
OLD( I  I )=TEMP 
11*11-1 
IF(X)85t85»90 

85  STRMX=(-SrCOF*THETA) 
PRINTl»TIME.TBARtSTRMX 
PRINTltALAM(  1)  #ALAM(2)  •GRB(l)  »GRB(2) 
IF(SENSE  SWITCH  3)91.89 

89  IF< (TEMP-PREV)-15.)92.92»3 

91  PUNCHltALPHA»TO.AKAP»COND 
PUNCH1.EPSB.DIV»STC0F»HR 
PUNCH l»ELSTBtALPHB» CON DB.DENSB 
PUNCH1»ELSTM.ALPHM»C0NDM»CP 
PUNCHl»ELSTtV»TIME.EPSM 
PUNCH1»DELT.HRSTD»TEMP 
D097I*1»4 

97  PUNCH1.QUK  I  )  tT(  I  ) 
PUNCH250»MtLL 
IF(SENSE    SWITCHD98.999 

98  D094I*1»8 

94  PUNCH1»FLUX(  I) 
D096I»1»11 

96  PUNCH1.0LD( I ) 

IF(SENSE  SWITCH  4)89»999 

92  DELT*2«0*DELT 

IF( (STDY-TEMP)-10«>95.3.3 

95  M»0 
LL«0 
TEMP-TO 
OUT(l)«TO 
TIME«0. 
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DELTaO.Ol 

GO  TO  2 
90  X=X-(EL/DIV) 

GO  TO  61 
100  HR2=EPS*SIGMA*( (TEMP**4)-( T0**4) )/ (TEMP-TO) 

IF(1-MM)115.115.102 
102  IF(2-KK)llu. 110.105 

105  IF( l-LL)lu6.1u6.158 

106  HRA=HR 
DEL1=HR-HR2 
DELSV=DEL1 
IF(DELl)107»120tl09 

107  DEL1=(-DEL1> 

109  IF(DEL1/HR-5.E-O5)120,108tl08 

108  KK  =  2 

GO  TO  160 

110  HRB*HR 
DEL2=HR-HR2 
DELL=DELSV+DEL1 
IF(DELL)113. 112.113 

112  DEL1=0ELSV 
IF(HRA-HRB)113.120tll3 

113  SL0PE=(DEL1-DEL2)/(HRA-HRB) 
IF(SLOPE)114.120.114 

114  B=DEL1-SL0PE*HRA 
HRSAV=(-B/SLOPE) 
KKK»1 

KK»1 

GO  TO  160 

115  IF(HRSTD-HR2) 116.116.117 

116  HR«(HR+HRSTl>)/2. 
GO  TO  118 

117  HR-HR2 

118  MM*0 

GO  TO  4 

120  PRINT1.HR 
IF(2-M)121.119.80 

119  OUT(2)*TEMP 

121  IF(3-M)123.122.80 

122  OUT(3)*TEMP 

123  IF(4-M) 124.124.80 

124  0UT(4)»TEMP 
GO    TO    80 

158    HR-HR2 

160  TBAR2*0. 
IFU-LD161.161.164 

161  D0162J«1»2 

TRM5«6RB( J)*SIN(VALUE( J) ) /VALUE (J) 

162  TBAR2*TBAR2+TRM5 
164  TBAR1*0« 

D0163I«1.8 

TRM11*FLUX< I )*( l«/HR-l./(COND*EL*U(I )**2) ) 

TRM1«(1.-EXP(-U( I >*EU )*TRM11 
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TRM2=FLUX( 1)*EXP(-U( I )  *EL) /(COND*U(i ) ) 
TRM3=FLUX( 1 ) *EL/ ( 2 .*COND ) 
163  TBAR1=TBAR1+TRM1+TRM2+TRM3 
THETA=1BAR1-TBAR2 
TBAR=THETA+TO 
IF(l-LL)172tl72»168 

168  rBVAL=TBSAV-TBAR 
IF<TBVAL)169»180»171 

169  TBVAL*(-TBVAL) 

171  IF( (TBVAL/TBAR)-l.E-06)180tl80>172 

172  TBSAV-TBAR 
COND=CONDB+CONDM* ( TBAR-460 • 0  ) 
ALPHA=ALPHB+ALPHM*( TBAR-460. ) 
DENS*DENSB/< 1.0+3.*ALPHA*( TBAR-50O.U) ) 
AKAP=COND/  (  OENS*CP  ) 
EPS=EPSB+EPSM*( TBAR-500. ) 

IF(TBAR-5  00.0)173»175»175 

173  ELST=ELSTB 
GO  TO  174 

17  5  ELST«ELSTB+ELSTM*(TBAR-560.0) 

174  STCOF=ELST*ALPHA/(1.0-V) 
IF(l-LLU130tl30»181 

181     IF(l-K.KM177tl77»176 

176  IF(2-KK)178»178tl79 

177  HR=HRSAV 
KKK  =  0 

GO  TO  179 

178  HR=HR2 

179  IF(1-LL)4»4»78 

180  PRINT1»HR 
HRSTD=HR 
GO  TO  70 

300  MM«1 

IF(2-M)310t31u»305 
305  T(2)*TIME 

M*M+1 

MM»0 

GO  TO  4 
310  TEMP«OUT(  1)+(TIME-T(  1)  )  *  (OUT  I  2  )  -OUT  (  1))/(T(2)-T(D) 

IF(3-M)320,320t313 
313  T(3)*TIME 

M»M+1 

GO  TO  100 
320  TMP21«(OUT(3)-OUT(2)  )/<T(3)-T(2>  )- (OUT ( 2 ) -OUT (  1))/(T(2)-T(1)) 

TEMP»TEMP+TMP21*(TIME-T(1>  )*(TIME-T12)  )/(T(3)-T(l)) 

IF(4-M)33O»330t323 
323  T(4)*TIME 

M«M+1 

GO  TO  100 
330  TMP31«(GUT(4)-QUT(3) )/(T(4)-T(3) ) - (OUT ( 3 ) -OUT ( 2 ) 1/ ( T( 3) -T( 2 ) I 

TMP32*(1./(T(4)-T(2) ))*TMP31 

TMP33«(1./(T(3)-T(2) ) )»TMP21 

TMP34»(TIME-T(1)  )*(TIME-T(2)  )*(TIME-T(3)  )/(T(4)-T(l)) 
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TEMP*TEMP+TMP34*(TMP32-TMP33  ) 

D0335I»1»3 

J*I+1 

OUT( I )«OUT( J) 
335  T(  I  )»T< J) 

T(4)*TIME 

0UT(4)»TEMP 

GO  TO  100 
400  PRINT1»STDY,X 

IF(X)402»402»405 
402  LL»1 

GO  TO  300 
405  X*X-(EL/DIV) 

GO  TO  70 

END 
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APPENDIX  E 


Description  and  Explanation  of  the  IBM  1620  FORTRAN 
Program  to  Calculate  the  Temperature  Distribution  and 
Thermal  Stresses  in  a  Plate  with  Insulation  at  One  Face 
and  Radiative  Heat  Loss  at  the  Other  Face,  when  Subjected 
to  a  Time  Dependent  Gamma  Photon  Beam  at  the  First  Face 


The  presence  of  the  time  dependent 
heat  generation  term  in  this  system 
makes  it  expedient  to  resort  to  the  use 
of  numerical  techniques  in  its  solution. 
The  partial  differential  equation  de- 
scribing this  system  is 


(E-1) 


9*  '     A      ajat 

FIG.  19. 

In  finite  difference  form,    this  equation     TEMPERATURE      DISTRIBUTION 


is  written  as 


(B-2) 


-»■<+* 


->M 


p  «.*»*» 


-   J_ 


which  may  be  solved  for  the  (n  +  1 )st  time  increment  at  each  of  the  j  mesh 
points  across  the  thickness  of  the  plate  as  shown  in  Fig.   19.     The  value  of 
4-*  lies  between  zero  and  one  and  gives  an  optimum  result  when  chosen  as   (18) 


T- 


Z.         \2.0}At 


(3-3) 


Analysis  of  the  resulting  system  of  equations  yields  the  matrix  found  in 
Appendix  B. 


103 


This  matrix  can  be  solved  by  any  one  of  a  number  of  methods.  The  one 
chosen  was  the  Crout  reduction.  A  special  subroutine,  called  CRAM,  has  been 
written  to  solve  an  augmented  matrix  by  the  Crout  reduction.   The  subroutine 
carries  20  significant  digits  in  its  solution  and  yields  answers  which  are 
quite  accurate  in  the  eight  significant  digits  reported  in  the  FORTRAN  out- 
put. The  use  of  this  extra  accuracy  was  deemed  necessary  when  an  ordinary 
Crout  reduction  was  found  to  have  round  off  errors  whose  magnitude  increased 
as  the  iteration  progresses.  It  was  found  that  after  two  or  three  time  in- 
crements the  answers  obtained  by  normal  methods  were  completely  meaningless. 

In  the  solution  of  the  temperature  distribution  at  each  time  increment, 
the  procedure  is  much  the  same  as  was  used  in  the  soltuion  to  Bq.(D-l).  A 
value  of  hr  is  first  assumed.   The  matrix  was  solved  for  the  outside  temper- 
ature and  hr  was  then  checked  against  the  radiation  boundary  condition.  The 
value  of  hr  found  from  this  boundary  condition  was  used  as  the  next  estimate 
and  was  again  checked.   Then  a  linear  approximation  was  used  as  the  first 
of  the  next  two  assumptions  for  hr.  When  the  assumed  value  came  within  the 
specified  proximity  of  the  boundary  condition,  the  iteration  was  considered 
complete. 

The  insertion  of  the  constants  into  the  matrix  first  requires  the  solu- 
tion of  Bq.(B-2)  to  find  the  gamma  photon  flux  at  the  time  of  interest. 
Since  many  of  the  matrix  coefficients  are  temperature  dependent,  it  is  nec- 
essary to  alter  them  after  each  trial  for  h  has  been  completed,  in  the  same 
manner  as  was  done  in  Appendix  D. 

When  the  true  value  of  hr  had  been  found,  the  free  plate  thermal  stress 
was  calculated  by  the  use  of  Simpson's  rule,  as  shown  in  Bq.(B-2l).  The 
temperature,  the  free  plate  thermal  stress  and  the  distance  from  the  origin 


1  (fl- 


are printed  as  output  after  the  value  for  hr  has  been  printed.  The  last 
printed  line  contains  the  time,  the  average  temperature  and  the  average  max- 
imum thermal  stress  for  a  perfectly  restrained  plate.  The  choice  as  to  when 
the  calculations  were  complete  was  completely  arbitrary,  depending  on  how 
closely  each  case  was  desired  to  be  investigated.  Consequently,  switch  2 
may  be  turned  ON  at  any  time  to  terminate  the  program  at  the  completion  of 
the  calculation  then  in  progress.  After  completion,  the  program  will  again 
start  at  the  beginning  with  a  new  set  of  input  data.  Switch  1  was  also  in- 
cluded in  the  program  so  that  the  time  increment  could  be  doubled  at  any  de- 
sired time  by  turning  the  switch  ON. 

The  entire  input  to  this  program  consists  of  the  punched  card  output 
from  the  program  described  in  Appendix  D.  A  single  card  is  put  at  the  front 
of  the  input  data  and  contains  the  absorption  coefficient  at  1  .05  Mev.  for 
the  material  being  examined.  The  size  of  tne  plate  being  considered  is  put 
on  a  data  card  and  follows  the  remainder  of  the  input  data. 

The  flow  chart  for  this  program  and  the  FORTRAN  source  deck  statements 
appear  on  the  following  pages. 
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106 


AL  =  AL/2 
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DIMENSION  SAV(  lOtll)  »FOR(  11)  »OLD(  11)  »FLUX(8)  tU(8) tT<4)  tOUT(4) 
1  FORMAT(E15.8»E15.8,E15.8,E15.8 ) 
READltU 

3  READl»ALPHA,TOtAKAP»COND 
READ1.EP5B»0IV.STC0F  »HR 
READ1»ELSTB»ALPHB»CQNDB»DENSB 

READl»ELSTM»ALPhM»CONDM»CP 
READl»ELSTtVtTIMEtEPSM 
READ1»DELT,HRSTD»TEMP 
D04I=lt8 

4  READ1»FLUX( I ) 
D05I=1»11 

5  READ1»0LD( I ) 
TAU=TIME*360Q. 
Z  =  0. 

EPS=EPSB 

SIGMA=1.7123E-09 
READ1»EL 
DELX=EL/DIV 
FLUXX=Q. 
D07I=1»8 

7  FLUXX=FLUXX+FLUX( I ) 

GO  TO  300 
12  TIME=TIME+DELT 

PREV=OLD( 11) 

Z=Z+DELT*3600. 

ZT=Z+TAU 

FFLUX=.1*( 1./UZ  +  10.)**(.2))-1./((ZT+1Q.)**1.2))) 

FFLUX=FFLUX+.08  7*( l./((Z+2.E+07)**l.2))-l./((ZT+2.E+07)**(.2))) 

FFLUX=FFLUX-.u025*(EXP(-Z/2040. ) -EXP ( -ZT / 2040. ) ) 

FFLUX=FLUXX*(FFLUX-.0013*(EXP(-Z/2  90000. ) -EXP ( -ZT/ 290000. ) ) ) 
15  BNDRY=( l.+HR*DELX/COND) 

PSI=Q.5-DELX**2/(12.*AKAP*DELT ) 

RAT I0=DELX**2/ ( DELT*AKAP ) 

RAPS1=RATI0+2.*PSI 

DIFF=1.-PSI 

FLC0N=U*FFLUX*DELX#*2/C0ND 

D030I=1»10 

D030J=1»11 
30    SAVU  »J)=0. 

DO20I=2»10 

SAV( I »11)=DIFF*(0LD( 1  +  1 )-2.*0LD( I )+OLD( 1-1) ) +RAT 10* ( OLD ( I )-T0) 

H-I-l 
2  0    SAV( I tll)=SAV(I»ll)+FLCON*EXP(-H*U*DELX) 

S AVI  1  til) =DIFF*(2.*OLD(2)-2.*OLD( D  ) +RAT I  0* ( OLD ( 1)-T0)+FLC0N 

M  =  l 

SAV(1»1 )=RAPSI 

SAV(lt2)=-2.*P$I 

D040I=2t9 

SAV( I »I-1)=-PSI 

SAVU  tl )=RAPSI 
40    SAVd  tI+l)=-PSI 

SAV(10»9)=-PSI 

5AV(lO»10)=RATIO+(2.-l./BNDRY)*PSI 
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FOR(l)=CRAM(  10.  ) 
80    FQR(ll)=FOR( 10) /BNDRY 

TEMP  =  FOR(  1D+TO 
100    HR2=EPS*SIGMA*(  (TEMP**4)-( T0**4) )/ I  TEMP-TO) 

IF(l-MM)115tll5»102 
102    IF(2-KK)110tllO»106 

106  HRA=HR 
DEL1=HR-HR2 
DELSV=DEL1 

IF (DELI )107tl20.l0S 

107  DEL1=(-DEL1) 

10  9  IF(DELl/HR-l.E-05)12G,108tl08 

108  HR=HR2 
KK  =  2 

GO  TO  200 
110  HRB=HR 

DEL2=HR-HR2 

DELL=DELSV+DEL1 

IF(DELL)113»112tll3 

112  DEL1=DELSV 
IF(HRA-HRB) 113tl20tll3 

113  SL0PE=(DEL1-DEL2)/(HRA-HR3) 
IF(SL0PE)114,120tll4 

114  B=DEL1-SL0PE*HRA 
HR=(-B/SL0PE) 
KK=1 

GO  TO  200 

115  HR=HR2 
MM  =  0 

GO  TO  15 
120  PRINT1»HR 

M«0 

GO  TO  200 
160  TBAR=STR2+T0 

COND=CQNDB+CONDM*(TBAR-460.0) 

ALPHA=ALPHB+ALPHM*( TBAR-46Q.0 ) 

DENS=DENSB/( 1.0+3 .*ALPHA* ( TBAR- 500.0) ) 

AK>AP  =  COND/(DENS*CP  ) 

EPS=EPSB+EPSM*( TBAR-500.) 

IF(TBAR-500.0)173tl7  5tl7  5 

173  ELST=ELSTB 
GO  TO  174 

175  ELST=ELSTB-ELSTM*( TBAR-500. 0) 

174  STC0F=ELST*ALPHA/(1.-V) 
GO  TO  15 

190  DEL=0. 

D0192I=ltll 

X=DEL*DELX 

DEL=DEL+1.G 

OLDd  )=F0R(  I  )+T0 
192    PUNCHltOLD(I)»STRES(I)tX 

IF(2-N)182.181»193 

181  OUT(2)=TEMP 

182  IF(3-N)184.183»193 
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183  OUT(3)=TEMP 

184  IF(4-N) 185tl85»193 

185  0UT(4)«TEMP 

STMAX=STR2*STCOF 

193  PUNCHl»TIMEf TBAR»STMAX 
IF(SENSE  SWITCH  2)3»250 

250  IFISENSE  SWITCH  1)251»12 

251  DELT=2«*DELT 
GO  TO  12 

200  STR2=F0R( 1 )+4.» ( FOR( 2 ) +FOR ( 4 ) +FOR ( 6 ) +  FOR { 8 ) +FOR ( 10) ) 

STR2=(STR2+2«*(FOR(3)+FOR( 5)+FOR( 7)+FOR(9) )+FOR( 11) )/30. 

IF(1-M)160, 160,202 
202  STR31=(FOR(ll)-FOR(l) )*EL/2. 

STR32=4.*( (9.*DELX-EL/2«)*FOR(2)+(DELX-EL/2«)*FOR( 10) ) 

STR33=2«*( (8.*DELX-EL/2. ) *FQR ( 3 ) + I 2.*DELX-EL/2. )*FOR(9) ) 

STR34=4.*( (7.*DELX-EL/2.)*FOR(4)+(3.*DELX-EL/2. )*FOR(8) ) 

STR35=2«*( (6.*DELX-EL/2.)*FQR(5)+(4.*DELX-EL/2«)*FOR(7) ) 

STR37=STR31+STR32+STR33+STR34+STR35 

R«0, 

DO205I=l»ll 

STRESl I )=STC0F*(STR2-F0Rl I ) +  .4* ( R*DELX-EL/2. ) *STR3  7/EL**2 ) 
205  R=R+1« 

GO  TO  190  . 
300  MM=1 

IF(1-N)310t3l0.305 
305  T(2)«TIME 

N*N+1 

MM«0 

GO  TO  15 
310  TEMP=OUT( 1)+(TIME-T( 1 ) ) * (OUT  I  2 ) -OUT ( 1))/(T(2)-T(1)) 

IF(2-N)320»320t313 
313  T(3)=TIME 

N  =  N+1 

GO  TO  100 
320  TMP213(0UT(3)-0UT(2)  )/ (T(3)-T(2)  )-(OUT  (2) -OUT  I  1))/(T(2)-T(D) 

TEMP  =  TEMP  +  TMP21*(TIME-T(1)  )*(TIME-T(2)  )/(T<3)-T(l)) 

IF(3-N)330»330»323 
323  T(4)=TIME 

N*N+1 

GO  TO  100 
330  TMP31-(OUT(4)-OUT(3) )/<T(4)-T(3) ) - ( OUT ( 3 ) -OUT ( 2 ) ) / ( T ( 3 ) -T ( 2 ) ) 

TMP32=(1«/(T(4)-T(2) ) )*TMP31 

TMP33=(1./(T(3)-T(2) ) )*TMP21 

TMP34=(TIME-T(1)  )*(TIME-T(2)  )*(TIME-T(3))/(T(4)-T(1)) 

TEMP=TEMP+TMP34*(TMP32-TMP33) 

•D0335I=1»3 

J*I+1 

OUT( I )«OUT( J) 
335  T( I )*T< J) 

T(4)«TIME 

0UT(4)«TEMP 

GO  TO  100 

END 
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APPENDIX  F 

An  Explanation  and  Description  of  the  IBM  1620  FORTRAN 
Program  to  Equate  the  Steady  State  Maximum  Thermal  Stress 
and  Find  the  Temperature  Distribution  in  the  Outer  Plate 

The  consideration  used  to  equate  the  maximum  thermal  stress  in  the  two 
plates  was  an  equality  of  the  average  temperatures  within  certain  tolerance 
limits,  which  were  arbitrarily  chosen.  The  equation  for  the  average  steady 
state  temperature  in  the  first  plate  is 

In  the  outer  plate  the  average  steady  state  temperature  is  given  by 

^=^~^t~"   ^^^^^  (F"2) 

where 


(F-3) 


'      _      KK  &  (l  -  *~u)  + 1,  (h^S-  K) 

and  also  where 

1,-1.^  (WO 

The  maximum  thermal  stress  for  the  completely  restrained  plate  is  given 
by 
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ow--^§-(8)  (F-5) 

and  in  the  present  case  this  has  been  altered  slightly  to  give  the  average 
maximum  thermal  stress.  This  is 

£.,=  -- £§-(8)  (F-6) 

so  that  a  single  uniform  stress  is  present  at  the  ends  of  the  plate.  It  is 
easily  seen  that  the  maximum  stress  in  each  plate  is  equal  when  the  average 
temperature  in  the  plates  are  equal.  The  two  temperatures  were  considered 
equal  when  the  quantity  AT/T  is  less  than  10"^. 

The  steady  state  temperatures  are  a  function  of  the  thichness  of  each 
of  the  two  plates  for  some  fixed  total  thickness  of  the  two  plates.  Hence, 
it  is  necessary  to  iterate  for  the  correct  thickness.  This  is  done  by  as- 
suming some  initial  percentage  of  the  total  thickness  is  in  the  first  plate. 
If  the  average  temperature  in  the  first  plate  is  higher  than  that  of  the  sec- 
ond plate,  its  thickness  is  decreased  by  some  fixed  increment  until  the  av- 
erage temperature  in  the  first  plate  is  less  than  that  of  the  second. 

At  this  point,  the  increment  is  halved  and  the  thickness  of  the  plate 
is  increased  until  the  temperature  is  again  greater  in  the  first  plate.  The 
increment  is  again  halved  and  the  iteration  then  proceeds  in  the  opposite 
direction.  This  procedure  continues  until  the  limits  have  been  reached. 
If  the  initial  estimate  yields  a  lower  average  temperature  in  the  first  plate, 
the  iteration  proceeds  in  the  opposite  order  from  that  described  above. 

When  the  iteration  has  been  completed,  the  temperature  distribution  in 
the  second  plate  is  calculated.  It  was  not  thought  necessary  to  calculate 
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the  distribution  in  the  first  plate  since  this  is  essentially  the  same  as 
the  case  of  the  single  plate  previously  considered.  The  temperature  distrib- 
ution in  the  second  plate  is  given  by 

where  R-|  ,  R?,  and  Ro  are  the  same  as  defined  above. 

The  input  data  cards  to  this  program  consist  of  the  following: 
1  .  The  flux  and  absorption  coefficient  for  each  of  the  eight  energy- 
groups  employed.   (8  cards) 

2.  The  intercept  and  slope  of  the  emissivity  with  regard  to  its  tem- 
perature dependency,  the  value  of  the  quantity  ocE/(1  -~P)   and  the 
base  temperature.   (1  card) 

3.  The  temperature  dependent  intercepts  of  the  modulus  of  elasticity, 
the  coefficient  of  thermal  expansion  and  the  conductivity,  and  the 
density  of  the  material  at  40°F.   (1  card) 

k.     The  temperature  dependent  slopes  of  the  modulus  of  elasticity,  the 
coefficient  of  thermal  expansion  and  the  conductivity,  and  the 
value  of  the  heat  capacity.   (1  card) 

5.  The  initial  values  of  the  modulus  of  elasticity,  the  coefficient 
of  thermal  expansion  and  the  conductivity,  and  Poisson's  ratio. 
(1  card) 

6.  The  heat  transfer  coefficient  between  the  plates,  and  an  initial 
estimate  of  the  equivalent  gray  body  convective  radiation  coeffi- 
cient, the  total  thickness  of  the  two  plates  and  the  steady  state 
temperature  of  a  single  plate  of  the  same  total  thickness.   (1  card) 

The  output  data,  in  printed  form,  consists  of  the  following: 
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1  .  The  percentage  of  the  total  thickness  contained  in  the  second  plate, 
the  final  equivalent  convective  radiation  coefficient  and  the  con- 
vection coefficient  between  the  two  plates. 

2.  The  average  temperature  of  the  first  and  second  plates  and  the  aver- 
age maximum  thermal  stress  of  both  plates. 

3.  The  temperature  profile  of  the  second  plate. 

It  is  also  possible  to  follow  the  path  of  the  difference  between  the 
average  temperatures  in  the  first  and  second  plates  by  turning  switch  2  ON. 
At  the- completion  of  each  calculation  for  a  given  choice  of  thicknesses  for 
each  of  the  two  plates,  this  difference  will  be  printed. 

The  flow  chart  for  this  program  and  the  FORTRAN  source  deck  statements 
appear  on  the  following  pages. 
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1  FORMAT ( El  5, 8 • El  5.8 tE 15.8 tE 15.8) 

DIMENSION  FLUXt 16) tU(8) 

D015I=1»8 
15  READ1»FLUX< I ) »U( I ) 

READl»EPSB»EPSM»STCOFtTO 

READl»ELSTB.ALPHBfCONDBtDENSB 

READltELSTMtALPHMtC.ONDMtCP 

READl.ELSTt ALPHA tCONDtV 
17  READl»H2tHRtELTOTtSTDY 

EL=.5*ELTOT 

DELTL=ELTOT/20. 

SIGMA=1.7123E-09 

L»0 

LL  =  0 

M«0 

KK=1 

TBSAV=0. 

TBR22*0. 
20  D025I«lt8 

K«I+8 
25  FLUXU)=FLUX( I )*EXP(-U< I )*EL) 

164  TBAR1=0. 
D0163I=1.8 

TRM1 INFLUX ( I )*( l./H2-l./(C0ND*EL*Ul I )**2) ) 
TRM1=(1.-EXP(-U( I )*EL) )*TRM11 
TRM2=FLUX( I )*EXP(-U( I ) *EL ) / ( C0ND*U ( 1 ) ) 
TRM3=FLUX( I ) *EL/ ( 2.*C0ND ) 
163  TBAR1=TBAR1+TRM1+TRM2+TRM3 
TBAR=TBAR1+T0 

168  TBVAL=TBSAV-TBAR 
IF(TBVAL)169.180»171 

169  TBVAL=(-TBVAL) 

171  IF(TBVAL/TBAR-l.E-06)180tl80tl7  2 

172  TBSAV=TBAR 

165  COND=CONDB+CONDM*( TBAR-460.0) 
IF(COND) 166tl66.167 

166  COND=0.1 

167  ALPHA=ALPHB+ALPHM*( TBAR-460. ) 
DENS=DENSB/( 1. 0+3.*ALPHA* ( TBAR-500.0) ) 
AKAP*COND/(DE/^S*CP  ) 
EPS»EPSB+EPSM*(TBAR-500. ) 
IF(TBAR-500.0)173»17  5tl75 

173  ELST=ELSTB 
GO  TO  174 

17  5  ELST=ELSTB+ELSTM*(TBAR-560.0) 

174  STCOF=ELST*ALPHA/(1.0-V) 
IF(l-LL)176tl76»164 

176  lF(2-KK)178.178t400 
178  HR-HR2 

GO  TO  400 
180  STRS1=(-STC0F*TBAR1) 

EL«ELTOT-EL 

X*EL 
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LL  =  1 

GO  TO  400 
100  HR2=EPS*SIGMA*( (TEMP**4)-l TQ**4) )/ (TEMP-TO) 
102  IF(2-KK)110»110.106 

106  HRA=HR 
DEL1=HR-HR2 
TBR21=TBAR2 
TBR2=TBR21-TBR22 
1F(TBR2)107.120»109 

107  TBR2=(-TBR2) 

109  IF(TBR2-5.E-02)120f 108»108 

108  KK.=2 

GO  TO  165 

110  HRB=HR 
DEL2=HR-HR2 
TBR22=TBAR2 
IF(HRA-HRB) 113» 12  0*113 

113  SL0PE=(DEL1-DEL2)/(HRA-HRB) 
IF(SLOPE)114»120tll4 

114  B»DEL1-SL0PE*HRA 
HR*(-B/SLOPE) 
KK«1 

GO  TO  165 
120  LL»0 

GO  TO  450 
400  TEMP»0. 

TBAR2*0. 

DO405I»lt8 

K=I+8 

UL=(-U(I)*EL) 

Rl=(-FLUX(i<;)/(COND*U(  I  )  )  ) 

R2N=(-R1*HR*H2*( l.-EXP(UL) ) +FLUX { K ) * ( HR+H2*EXP ( UL ) ) ) 

R2=(-R2N/(COND*(HR+H2)+HR*H2*EL) ) 

R3=(-R2*<COND/HR+EL) ) +FLUX ( K ) *EXP I UL ) *  I l./(COND*U( I ) )-l./HR) 

TEMP=TEMP+R1*EXP(-U( I )*X)+R2*X+R3 

405  TBAR2=TBAR2+R1*(EXP(UL)-1. ) /UL+R2*EL/2.+R3 
TEMP=TEMP+TO 
IF(TBARl-TBAR2-TBARl/2.)40  4»464.464 

404  IF(TBAR2+T0-STDY>406»406»407 

406  TBAR=TBAR2+T0 
GO  TO  408 

407  TBAR=STDY 

408  IF(1-L)505.505»100 

450  TBR=TBAR1-TBAR2 
IF(SENSE  SWITCH2)449t448 

449  PRINT1»TBR 
448  TBSV2«TBR 

IF(TBR)451»500»452 

451  TBR«(-TBR) 

452  IF(TBR/TBAR2-l.E-04)500t500.453 

453  IF(TBSV2)455t500t465 

455  IF(1-M)475.456»456 

456  EL*ELTOT-EL+DELTL 
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464 


465 
466 


475 

500 


505 
506 

507 


M=l 

GO 

LL* 

EL* 

GO 

IF( 

EL  = 
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ABSTRACT 

In  considering  the  structural  aspects  of  nuclear  rockets,  it  is  possible 
to  construct  idealized  models  of  structures  which  might  be  found  in  the  vicin- 
ity of  the  reactor.  Such  a  model  is  an  infinite  flat  plate  one  face  of  which 
is  subjected  to  a  constant  gamma  photon  flux.  The  temperatures  and  thermal 
stresses  of  both  a  steel  and  an  aluminum  plate  are  investigated.   The  tempera- 
tures of  steel  plates  as  thick  as  one  inch  are  well  below  the  melting  point 
for  a  flux  intensity  of  9^-500  Btu/hr,  ft2  when  radiation  heat  loss  is  the  only 
mode  of  heat  transfer  available.   For  aluminum,  this  energy  flux  must  be  re- 
duced by  a  factor  of  ten  to  keep  it  below  the  melting  point. 

In  steel,  the  maximum  thermal  stress  for  a  plate  held  in  perfect  re- 
straint exceeds  the  yield  stress  after  being  subjected  to  the  gamma  heating 
for  a  few  minutes.  The  aluminum  must  be  heated  for  longer  than  one  hour  be- 
fore the  maximum  thermal  stress  exceeds  the  yield  stress.   For  a  plate  re- 
strained in  bending,  the  thermal  stresses  in  both  aluminum  and  steel  are  be- 
low the  yield  stress.   The  free  plate  stresses  are  small  enough  to  be  consid- 
ered insignificant. 

After  shutdown  of  the  reactor,  the  plate  is  exposed  to  a  time  dependent 
gamma  flux.   Although  its  magnitude  is  quite  significant,  the  plate  initially 
cools  at  a  rate  almost  the  same  as  if  no  additional  heating  were  present. 
The  cooling  begins  so  rapidly  that  the  temperature  appears  to  be  discontinu- 
ous with  time. 

Additional  cooling  can  be  provided  for  the  plate  by  constructing  a  sys- 
tem which  enables  the  other  face  of  the  plate  to  lose  heat  by  convection. 
Two  parallel  plates  were  considered,  where  the  total  thickness  of  both  plates 
was  the  same  as  that  of  the  single  aforementioned  plate.   Then  the  steady 


state  temperatures  and  thermal  stresses  were  compared  with  the  single  plate 
case  as  the  magnitude  of  the  convection  coefficient  was  changed.  The  base 
temperature  to  which  heat  was  lost  by  both  radiation  and  convection  was  al- 
ways assumed  to  be  500°R. 

The  temperature  in  both  steel  and  aluminum  were  found  to  be  lower  than 
the  single  plate  case.   In  aluminum,  this  difference  was  most  pronounced 
since  it  has  a  low  emissivity  and  does  not  readily  lose  heat  by  radiation 
alone.   The  thermal  stress  for  the  case  of  a  plate  restrained  in  bending  ex- 
hibited a  minimum  at  the  time  when  the  heat  lost  from  one  face  of  the  plate 
by  radiation  was  equal  to  the  heat  lost  from  its  other  face  by  convection. 
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